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Is the Universe logotropic? 
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We consider the possibility that the universe is made of a single dark fluid described by a logotropic 
equation of state P = A\n(p/p,), where p is the rest-mass density, p* is a reference density, and A is 
the logotropic temperature. The energy density e is the sum of two terms: a rest-mass energy term 
pc? that mimics dark matter and an internal energy term u(p) = —P(p) — A that mimics dark energy. 

This decomposition leads to a natural, and physical, unification of dark matter and dark energy, and 
elucidates their mysterious nature. In the early universe, the dark fluid behaves as pressureless dark 
matter (P ~ 0, e oc a -3 ) and, in the late universe, it behaves as dark energy ( P ~ —e, t oc In a). The 
logotropic model depends on a single parameter B = A/pA<? (dimensionless logotropic temperature) 
where pA = 6.72 x 10 -24 gm -3 is the cosmological density. For B = 0, we recover the ACDM model 
with a different justification. For B > 0, we can describe deviations from the ACDM model. Using 
cosmological constraints, we find that 0 < B < 0.09425. We consider the possibility that dark matter 
halos are described by the same logotropic equation of state. When B > 0, pressure gradients prevent 
gravitational collapse and provide halo density cores instead of cuspy density profiles, in agreement 
with the observations. The universal rotation curve of logotropic dark matter halos is consistent 
with the observational Burkert profile up to the halo radius. It decreases as r~ 1 at large distances, 
similarly to the profile of dark matter halos close to the core radius [Burkert, arXiv:1501.06604 . 
Interestingly, if we assume that all the dark matter halos have the same logotropic temperature 
B, we find that their surface density E = porh is constant. This result is in agreement with the 
observations [Donato et al, MNRAS 397 , 1169 (2009)] where it is found that Eo = 141 Alg/pc 2 for 
dark matter halos differing by several orders of magnitude in size. Using this observational result, 
we obtain B = 3.53 x 10 -3 . Then, we show that the mass enclosed within a sphere of fixed radius 
r u = 300pc has the same value M 300 = 1-93 x 10‘ M@ for all the dwarf halos, in agreement with 
the observations [Strigari et al, Nature 454, 1096 (2008)]. Finally, assuming that p* = pp , where 
p P = 5.16 x 10" gm -3 is the Planck density, we predict B = 3.53 x 10 —3 , in perfect agreement with 
the value obtained from the observations. We approximately have B ~ l/ln(pp/pA) — l/[1231n(10)[ 
where 123 is the famous number occurring in the ratio pp/pA ~ 10 123 between the Planck density 
and the cosmological density. This value of B is sufficiently low to satisfy the cosmological bound 
0 < B < 0.09425 and sufficiently large to differ from CDM ( B = 0) and avoid density cusps in dark 
matter halos. It leads to a Jeans length at the beginning of the matter era of the order of A j = 40.4 pc 
which is consistent with the minimum size of dark matter halos observed in the universe. Therefore, 
a logotropic equation of state is a good candidate to account both for galactic and cosmological 
observations. This may be a hint that dark matter and dark energy are the manifestation of a single 
dark fluid. If we assume that the dark fluid is made of a self-interacting scalar field, representing 
for example Bose-Einstein condensates, we find that the logotropic equation of state arises from the 
Gross-Pitaevskii equation with an inverted quadratic potential, or from the Klein-Gordon equation 
with a logarithmic potential. We also relate the logotropic equation of state to Tsallis generalized 
thermodynamics and to the Cardassian model (motivated by the existence of extra-dimensions). 

PACS numbers: 95.30.Sf, 95.35.+d, 95.36.+X, 98.62.Gq, 98.80.-k 


I. INTRODUCTION 


The cosmological constant has an interesting history. 
Its pre-history dates back to 1896, when von Seeliger [T] 
and Neumann [2] introduced an attenuation factor in the 
Newtonian potential, $ = —Gme~ Xr /r, in order to re¬ 
solve the non-convergence of the gravitational force in 
Newtonian cosmology for an infinite homogeneous distri¬ 
bution of matter [3j. In 1917, Einstein [4] introduced a 
cosmological constant A in the equations of general rel¬ 
ativity in order to obtain an infinite homogeneous static 
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universe. 1 After the discovery of the expansion of the 
universe, predicted theoretically by Friedmann uni nu 


1 As a preamble of his paper, in order to motivate the cosmological 
constant, Einstein [4] considered the Newtonian approximation 
and replaced the Poisson equation by an equation of the form 
A<f> — (A/c 2 )4> = 47r Gp, so as to obtain a homogeneous static 
solution 4> = —AnGpc 2 /A. In Einstein’s belief, the Newtonian 
effect of the cosmological constant was to shield the gravitational 
interaction on a distance c/ vA, similarly to the proposal of Seel¬ 
iger and Neumann. This is actually incorrect. Lemaitre [5] was 
the first to understand that the cosmological constant can be 
interpreted as a force of “cosmic repulsion” (see also [6lf9]i. In 
the Newtonian context, the modified Poisson equation including 
the cosmological constant is A4> = AnGp — A, leading to the 
homogeneous static solution p = A/AnG. 
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and Lemaitre mm, and ascertained by the observa¬ 
tions of Hubble ! M] , Einstein considered the cosmologi¬ 
cal constant as “the biggest blunder of his life” he and 
banished it. In 1932, Einstein and de Sitter m pub¬ 
lished a short paper describing the expansion of a pres¬ 
sureless universe without cosmological constant. This 
is the so-called Einstein-de Sitter (EdS) model. In this 
model, the universe undergoes a decelerating expansion. 
The discovery of the acceleration of the expansion of the 
universe in recent years nzj revived the interest in the 
cosmological constant since a positive cosmological con¬ 
stant can precisely account for such an acceleration. * 2 At 
present, the standard model of cosmology is based on a 
pressureless dark matter fluid (P = 0) with a positive 
cosmological constant (A > 0). This is the cold dark 
matter (ACDM) model. From the observations, the cos¬ 
mological constant has the value A = 1.13 x 10 -35 s -2 . 
The cosmological constant A is equivalent to a constant 
energy density = Ac 2 /8nG = 6.04 x 10~ 7 gm -1 s -2 , 
called dark energy, associated with an equation of state 
P = — e involving a negative pressure. Some authors 
EM] have proposed to interpret the cosmological con¬ 
stant in terms of the energy of the vacuum. However, 
when one tries to explain the cosmological constant in 
relation to the vacuum energy, one is confronted to the 
so-called “cosmological constant problem” [5|; [55] be¬ 
cause quantum field theory predicts that the vacuum en¬ 
ergy density should be of the order of the Planck density 
Pp = c 5 * /HG 2 = 5.16 x 10" gm -3 which is 123 orders of 
magnitude larger than the value of the cosmological den¬ 
sity pa = ca/c 2 = 6.72 x 10 - 24 gm -3 deduced from the 
observations. To circumvent this problem, some authors 
have proposed to abandon the cosmological constant A 
and to explain the acceleration of the universe in terms of 
a dark energy with a time-varying density associated to 
a scalar field called “quintessence” [26]. This mysterious 
concept of “dark energy” adds to the other mysterious 
concept of “dark matter” necessary to account for the 
missing mass of the galaxies inferred from the virial theo¬ 
rem m and to explain their flat rotation curves [28j |29| . 
As an alternative to quintessence, Kamenshchik et al. 
[50] proposed a heuristic unification of dark matter and 
dark energy in terms of an exotic fluid with an equation 
of state P = — A/e , called the Chaplygin gas. This equa¬ 
tion of state provides a model of universe that behaves 
as a pressureless fluid (dark matter) at early times, and 
as a fluid with a constant energy density (dark energy) 
at late times, yielding an exponential acceleration simi¬ 
lar to the effect of the cosmological constant. However, 
in the intermediate regime of interest, this model does 
not give a good agreement with the observations mm ■ 


2 The effect of the cosmological constant in an empty universe 
(p = 0) was first considered by de Sitter |18l I19| in the context 
of a static universe. It was understood later by Lemaitre M\ 
that de Sitter’s solution actually describes a universe expanding 
exponentially rapidly with time. 


Therefore, generalizations of the Chaplygin gas have been 
considered [33]. We note that an equation of state corre¬ 
sponding to a constant negative pressure P = —e a also 
provides a heuristic unification of dark matter and dark 
energy, and is simpler [551-1571. Furthermore, this equa¬ 
tion of state is consistent with the observations since it 
gives the same results as the ACDM model although its 
physical justification is different. Therefore, we can re¬ 
alize a simple unification of dark matter and dark en¬ 
ergy by assuming that the universe is filled with a single 
dark fluid with a constant pressure. However, this model 
has two drawbacks. Since it has no free parameter, it 
cannot account for (small) deviations from the ACDM 
model. On the other hand, a constant equation of state 
cannot describe dark matter halos. Indeed, since there 
is no pressure gradient (VP = 0 ), this model behaves 
similarly to the CDM model where P = 0, and gener¬ 
ates cuspy density profiles at the center of the halos [38]. 
This is problematic because observations of dark matter 
halos favor density cores instead of cusps Emu]. Cores 
can be naturally explained by assuming that dark mat¬ 
ter halos are collisional and that they are described by 
an equation of state P{p) that strictly increases with the 
density so that pressure gradients can counterbalance the 
gravitational attraction. Therefore, if we want to realize 
a successful unification of dark matter and dark energy, 
the idea is to consider a dark fluid with an equation of 
state P{p) increasing slowly with the density. 

A natural idea is to consider a polytropic equation of 
state where the pressure is a power of the density with 
a small index 7 ~ 0. There are two types of polytropic 
equations of state in general relativity, as first realized 
by Tooper HU [42] in the context of general relativistic 
stars. 3 Model I corresponds to an equation of state of 
the form P = Ke 1 with 7 = 1 + 1/n, where the pressure 
is a power of the energy density. This model has been 
studied in detail in [3411571 in a cosmological context. It 
can be viewed as a generalization of the Chaplygin gas. 
In this model, when K < 0 and 7 < 1 (i.e. n < 0) 
[55], the dark fluid behaves as a pressureless fluid (dark 
matter) in the early universe and as a fluid with a con¬ 
stant energy density (dark energy) in the late universe. 4 
However, this distinction is only asymptotic and, in the 


! The cosmological implications of these two types of polytropic 

models have been recently discussed by the author in [431 HU- 
Similar arguments have been developed independently in m 

without reference to the work of Tooper. 

4 When K < 0 and 7 > 1 (i.e. n > 0) m, this polytropic model 

describes a phase of early inflation followed by a decelerating ex¬ 
pansion (modeling, e.g., radiation). In this model, the tempera¬ 

ture is initially very low, increases exponentially rapidly during 
the inflation up to the Planck temperature (GUT scale), then 

decreases algebraically. By contrast, in the usual inflationary 
scenario 1 461149 1 , the temperature decreases during the inflation 
and one must advocate a phase of re-heating (not very well un¬ 
derstood) to account for the observations. On the other hand, 
the polytropic model generates a value of the entropy as large 
as S/kg = 5.04 x 10 87 |34| . It is also possible to introduce a 
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intermediate regime, it is not possible to separate these 
two components unambiguously because they do not ap¬ 
pear additively in the energy density e unless 7 = 0 , 
which is equivalent to the ACDM model. Model II cor¬ 
responds to an equation of state of the form P = Kp 1 
with 7=1 + 1 /n, where the pressure is a power of the 
rest-mass density. In this model, the energy density e 
is the sum of two terms: a rest-mass energy pc 2 and an 
internal energy u(p) = P(p)/( 7 — 1). Since p ex a -3 , the 
rest mass term mimics dark matter. When K < 0 and 
7 < 1 (i.e. n < 0 ), the internal energy term mimics dark 
energy. This decomposition leads to a natural, and physi¬ 
cal, unification of dark matter and dark energy, and eluci¬ 
dates their mysterious nature. In order to account for the 
observations, these polytropic models must be relatively 
close to the ACDM model, equivalent to a constant neg¬ 
ative pressure. This implies that the polytropic index 7 
must be close to zero. For 7 = 0 (i.e. n = — 1), we recover 
the ACDM model exactly. For 7 different from zero, we 
can study deviations from the ACDM model. The poly¬ 
tropic models have been confronted to the observations 
in order to constrain the values of the admissible index 
7 or n. For Model I, it is found that n = —1.05+o '16 
[ST] . For Model II, it is found that —0.089 < 7 < 0 
(corresponding to — 1 < n < —0.92) [45]. Although these 
polytropic models are interesting in themselves, we antic¬ 
ipate a problem in order to justify them from a fundamen¬ 
tal theory. Actually, it is possible to justify a polytropic 
equation of state by assuming that the dark fluid is made 
of a self-interacting scalar field representing, for example, 
Bose-Einstein condensates (BECs). A polytropic equa¬ 
tion of state with 7 0 arises from the Gross-Pitaevskii 

(GP) equation with a power-law potential |^| 2 ( 7-1 ) or 
from the Klein-Gordon (KG) equation with a power-law 
potential |^>| 27 . However, it is hard to imagine how one 
can justify from fundamental arguments why the uni¬ 
verse is described by a polytropic equation of state with 
a “curious” index equal, for example, to 7 = —0.089. 
The most natural value of the polytropic index is 7 = 0 
corresponding to a constant pressure but, in that case, 
the GP and KG equations degenerate. Furthermore, a 
polytropic equation of state with 7 = 0 has no pressure 
gradient, which is problematic to account for the flat core 
of dark matter halos as previously mentioned. 

In this paper, we develop another idea that was 
sketched in Appendix B of 03}. We propose to describe 
the universe by a single dark fluid with a logotropic equa¬ 
tion of state of the form P = A\afp/p t ), where p is the 
rest-mass density, p* is a reference density, and A is the 
logotropic temperature. This equation of state was intro¬ 
duced phenomenologically in astrophysics by McLaugh¬ 
lin and Pudritz [52] to describe the internal structure and 


generalized polytropic model based on a quadratic equation of 
state that describes simultaneously the early inflation, the inter¬ 
mediate decelerating expansion, and the late acceleration of the 
universe [31- 


th.e average properties of molecular clouds and clumps. It 
was also studied by Chavanis and Sire [53] in the context 
of Tsallis generalized thermodynamics [5T where it was 
shown to correspond to a polytropic equation of state 
with 7 — > 0 and K —► 00 in such a way that A = 7 K is 
finite . 5 This limit is very relevant to the present cosmo¬ 
logical context because we precisely require an equation 
of state with a polytropic exponent close to 7 = 0. There¬ 
fore, the logotropic equation of state is a natural, and se¬ 
rious, candidate for the unification of dark matter and 
dark energy. Since the pressure must increase with the 
density, the logotropic temperature must be strictly pos¬ 
itive (A > 0) which is satisfying from a thermodynamical 
point of view. In the early universe, the dark fluid be¬ 
haves as pressureless dark matter (P ~ 0, e oc a~ 3 ) and, 
in the late universe, it behaves as dark energy (P ~ —e, 
e ex In a) . The logotropic model depends on a single pa¬ 
rameter B = A/pac 2 (dimensionless logotropic tempera¬ 
ture), where p\ is the cosmological density. For B = 0, 
we recover the ACDM model with a different justifica¬ 
tion. For B > 0, we can describe deviations from the 
ACDM model. Using cosmological constraints, we find 
that 0 < B < 0.09425. For B = 0.09425, the universe 
is normal for t < 20.7 Gyrs (the energy density decreases 
as the scale factor increases) and becomes phantom af¬ 
terwards (the energy density increases with the scale fac¬ 
tor). Actually, the logotropic model breaks down before 
entering in the phantom regime because the velocity of 
sound exceeds the velocity of light when t > 17.3 Gyrs. 
However, this moment is quite remote in the future. 

Since the logotropic equation of state is expected to 
provide a unification of dark energy and dark matter, 
it should also describe the structure of dark matter ha¬ 
los. In this respect, the logotropic equation of state has 
several nice properties. First of all, since the pressure in¬ 
creases with the density, pressure gradients prevent gravi¬ 
tational collapse and provide halo density cores instead of 
cuspy density profiles, in agreement with the observations 
mm- On the other hand, the universal rotation curve 
of logotropic dark matter halos is consistent with the ob¬ 
servational Burkert profile [39] up to the halo radius . 6 


5 In the framework of generalized thermodynamics, one can re¬ 
ally interpret A as a generalized temperature associated with a 
generalized entropy called Log-entropy [53] (see Appendix |C 5 c[ >. 

6 The logotropic rotation curve continues to increase at larger dis- 
tances while the Burkert rotation curve decreases. This feature 
comes from the fact that the logotropic density decreases as r -1 
while the Burkert density decreases as r ~ 3 . Interestingly, we 
note that the profile of dark matter halos close to the core radius 
decreases as r -1 m in agreement with the logotropic profile. 
On the other hand, some galaxies such as low surface bright¬ 
ness (LSB) galaxies have rotation curves that strictly increase 
with the distance. These galaxies are particularly well isolated 
and little affected by tidal effects. In other cases, tidal effects 
or complex physical processes (e.g. incomplete relaxation) have 
to be taken into account. This can affect the behavior of the 
density profile and of the rotation curve at large distances and 
explain how we pass from a r~ 1 behavior to a r~ 3 behavior. 
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Finally, and very interestingly, if we assume that all the 
dark matter halos have the same logotropic temperature 
B (a consequence of our model), we find that their surface 
density E = pofh is constant. This result is in agreement 
with the observations, where it is found that Eo has the 
same value Eo = 141M 0 /pc 2 for dark matter halos dif¬ 
fering by several orders of magnitude in size j56H58) . a 
result unexplained so far. Using this observational re¬ 
sult, we obtain B = 3.53 x 10 -3 . Then, we show that the 
mass enclosed within a sphere of fixed radius r u = 300 pc 
has the same value M 300 = 1.93 x 10' M 0 for all the 
dwarf halos, in agreement with the observations of Stri- 
gari et al. 1521 ■ This so far unexplained result is a direct 
consequence of the logotropic equation of state. Finally, 
assuming that p* = pp, where pp is the Planck density, 
we predict B = 3.53 x 10~ 3 , in perfect agreement with the 
value obtained from the observations. We approximately 
have B ~ 1/\n(pp/p\) — 1/[123ln(10)] where 123 is the 
famous number occurring in the ratio pp/pa ~ 10 123 be¬ 
tween the Planck density pp = 5.16 x 10" gm“ 3 and the 
cosmological density p\ = 6.72 x 10 -24 gm -3 . This value 
of B is sufficiently low to satisfy the cosmological bound 
0 < B < 0.09425 and sufficiently large to differ from 
CDM (B = 0) and avoid density cusps in dark matter 
halos. Therefore, a logotropic equation of state is a good 
candidate to account both for galactic and cosmological 
observations. This may be a hint that dark matter and 
dark energy are the manifestation of a single dark fluid. 

Finally, it may be less problematic to theoretically jus¬ 
tify the logotropic model than a polytropic model with 
an index 7 = —0.089 for example |45j. Indeed, the lo¬ 
gotropic equation of state has a universal functional form 
(logarithmic) while the polytropic model depends on an 
unknown index 7. As shown in detail below, the lo¬ 
gotropic model removes the degeneracy of the polytropic 
model when 7 = 0 . Therefore, there is some hope to jus¬ 
tify the logotropic model from a fundamental field the¬ 
ory. If we assume that the dark fluid is made of a self¬ 
interacting scalar field, representing for example BECs, 
we find that the logotropic equation of state arises from 
the GP equation with an inverted quadratic potential 
or from the KG equation with a logarithmic poten¬ 
tial In 10|. To our knowledge, these equations have not 
been considered previously. They are, of course, inter¬ 
esting in themselves and our approach strongly suggests 
that they may have important applications in cosmology. 

The paper is organized as follows. In Sec. |TTJ we re¬ 
view different cosmological models that have been pro¬ 
posed in the past. We start with well-known models 
(Einstein-de Sitter, ACDM, Chaplygin) then discuss less 
well-known ones (polytropic models of type I and II, Car- 
dassian model) that are important in our study. In Sec. 
m we introduce the logotropic model. In Sec. m we 
study the logotropic equation of state in a cosmological 
setting, as a unification of dark matter and dark energy, 
and obtain a cosmological bound 0 < B < 0.09425 for 
the normalized logotropic temperature B. In Sec. [Vj we 
apply the logotropic equation of state to the structure 


of dark matter halos and obtain B = 3.53 x 10 -3 from 
observations. In Sec. VI we predict B = 3.53 x 10 -3 
from theoretical considerations. In Appendix [Aj we pro¬ 
vide general results concerning the thermodynamics of a 
relativistic fluid. In Appendix [B] we provide a simple 
argument to justify the logotropic equation of state and 
mention its relation to a fifth force. In Appendix [Cj we 
develop a field theory for the dark fluid and determine 
the form of the GP and KG equations that generate a 
logotropic equation of state, as well as the corresponding 
generalized entropy. 

Remark: Secs. UH and m show the close connections 
between different cosmological models (many of these 
connections having not been discussed before) and ex¬ 
plain the physical reasons why we consider a logotropic 
equation of state. These sections are important to moti¬ 
vate our study. They also introduce our general formal¬ 
ism. However, the reader just interested by the results 
may skip these sections and go directly to Sec. [TV]for the 
application of the logotropic equation of state to cosmol¬ 
ogy and to Sec. [V] for the application of the logotropic 
equation of state to dark matter halos. 


II. COSMOLOGICAL MODELS 


A. The Friedmann equations 

We assume that the universe is homogeneous and 
isotropic, and contains a uniform perfect fluid of energy 
density e(t) and isotropic pressure P(t). The radius of 
curvature of the 3-dimensional space, or scale factor, is 
noted a(t) and the curvature of space is noted k. The 
universe is closed if k > 0 , flat if k = 0 , and open if 
k < 0. We assume that the universe is flat (k = 0) in 
agreement with the observations of the cosmic microwave 
background (CMB) [60]. In that case, the Einstein equa¬ 
tions can be written as m : 

— + 3-(e + P) = 0, (1) 

at a 


a 

a 


4t tG 

lk? 


(e + 3P) 


A 
3 ’ 


( 2 ) 


H 2 


fa\ 2 _8nG A 
\ a J 3 c 2 3 ’ 


(3) 


where we have introduced the Hubble parameter H = 
d/a and accounted for a possible non-zero cosmological 
constant A. 

Eqs. ([l])-([3]) are the well-known Friedmann equations 
describing a non-static universe. Among these three 
equations, only two are independent. The first equation 
can be viewed as an equation of continuity. For a given 
barotropic equation of state P = P(e), it determines the 
relation between the energy density e and the scale factor 
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a. Then, the evolution of the scale factor a(t) is given by 
Eq. (§. 

Introducing the equation of state parameter w = P/e, 
and assuming A = 0, we see from Eq. © that the uni¬ 
verse is decelerating if w > —1/3 (strong energy condi¬ 
tion) and accelerating if w < —1/3. On the other hand, 
according to Eq. the energy density decreases with 
the scale factor if w > —1 (null dominant energy condi¬ 
tion) and increases with the scale factor if it? < —1. The 
latter case corresponds to a “phantom” universe [62] . 


B. Einstein-de Sitter model 


We consider a universe made of pressureless matter 
(dust) with an equation of state 


P = 0, 


(4) 


and we assume that A = 0. This is the Einstein-de Sit¬ 
ter model 16]. The equation of continuity |T]) can be 
integrated into 




(5) 


where the subscript 0 refers to present-day values. The 
Friedmann equation (J3| writes 
2 


87rGeo 

3c 2 


(?)■ 


and it can be integrated into 

2/3 


a 

do 


= -Hnt 


e = 


67 t Gt 2 ’ 


( 6 ) 


(7) 


where H 0 = (SirGeo/Sc 2 ) 1 / 2 is the present value 
of the Hubble constant. From observations Hq = 
70.2kms _1 Mpc^ 1 = 2.275 10 -18 s -1 yielding eo = 
8.32 x 10 -7 gm _1 s“ 2 . In this model, the universe is al¬ 
ways decelerating and its age is 

2 


.EdS _ 

Lr\ — 


3Po 


( 8 ) 


Numerically, tg dS = 9.29 Gyrs. 


C. ACDM model 


D. Dark matter and dark energy 


As is clear from Eq. (|3]) , the cosmological constant is 
equivalent to a fluid with a constant energy density 7 


£A = Pac 2 


Ac 2 

87tG' 


( 10 ) 


According to Eq. ([!]), the equation of state leading to a 
constant energy density is P = —e. Instead of introduc¬ 
ing a cosmological constant A, we can consider that the 
universe is made of two fluids: a pressureless dark matter 
fluid with an equation of state P = 0 leading to a dark 
matter density e m = e mj o(ao/a) 3 , and a dark energy fluid 
with an equation of state 


P = -e (11) 

leading to a constant density identified with = p\c 2 = 
6.04 x 10~' gm -1 s -2 . In this context, is called the 
dark energy density, or the cosmological density. The 
total energy density of the universe (e = e m + £a) is 
therefore 8 


e = eo(~) + e A- ( 12 ) 

As a result, the Friedmann equation © can be written 
as 


H 

Ho 



+ ^, 0 ; 


(13) 


where we have introduced the present fractions of dark 
matter and dark energy fl m ,o = £m,o/eo and = 

ea/eo- By construction, o + 0 a, o = 1. From ob¬ 
servations = 0.274 and Ha,o = 0.726. Eq. (131 has 
a well-known analytical solution 


a 

do 


O \ 1 / 3 
W-M sinh 2/3 

^A ,0 / 



(14) 


£ _ Ha,o 

£q tanh 2 (§^/fiA,oH 0 t) 


(15) 


We still consider a universe made of pressureless matter 
(dust) with P = 0 but, in order to account for the accel¬ 
eration of the expansion of the universe, we now assume 
that A > 0. This is the ACDM model. The Friedmann 
equation (|3| now writes 

a\ 2 _ 87rGe 0 /a 0 \ 3 A 
a) ~ 3c 2 V a ) + 3' 1 j 

As discussed in the next section, this equation can be 
interpreted in terms of dark matter and dark energy, and 
it can be solved analytically. 


7 This energy was initially interpreted as the vacuum energy m- 
1231 . However, from particle physics and quantum field theory, 
the vacuum energy is expected to be of the order of the Planck 
density pp = 5.16 X 10 99 gm -3 while cosmological observations 
give pa = 6.72 x 10 —24 gm -3 . These densities differ by the huge 
ratio pp/pA ~ 10 123 . This is the so-called cosmological constant 
problem 1241125| . 

8 In this paper, we consider a sufficiently old universe so that radi¬ 
ation can be neglected. On the other hand, baryonic matter can 
be treated as a pressureless fluid (P& = 0) whose contribution 
adds to the contribution of dark matter. We can take it into ac¬ 
count in the previous expressions by considering that e m refers 
to dark matter and baryonic matter. 
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For early times, we recover the EdS solution (with mod¬ 
ified coefficients): 



e 4 
eo ~ 9 


(16) 


and for late times, we recover the de Sitter solution 


ao \4Ha,o/ 


e ~ e A . 


In the ACDM model, the age of the universe is 

2 


to = 


3tf 0 ^A.O 


sinh 


/ n A ,o x 1/2 ' 

\ C 


(17) 


(18) 


Numerically, f£ CDM = 13.8 Gyrs. 

In this model, the universe is decelerating at early 
times and accelerating at late times. The moment at 
which the universe starts accelerating (a > 0) takes place 
when e c = 3e A , he. a c /ao = (II mj o/2flA,o) 1 ^ 3 = 0.574. 9 
This corresponds to a time 

t c = - ^ sinh -1 ( ~^= ) . (19) 

\V2J 

Numerically, ^acdm _ 7 .18 Gyrs. The transition be¬ 
tween the dark matter era and the dark energy era takes 
place when e m = e A , hence £2 = 2e A , i.e. a^/a^ = 
(hl m) o/f! a.o) 1 ^ 3 = 0.723. This corresponds to a time 


t 2 = - . sinh 1 (1). 

3 H 0 y/n^ 

Numerically, fACDM = 9.61 Gyrs. 


( 20 ) 


E. The Chaplygin gas 


Kamenshchik et al. [501 have proposed to unify dark 
matter and dark energy by a single dark fluid described 
by the equation of state 

P=- A . (21) 

e 

This is called the Chaplygin gas because this equation of 
state was introduced by Chaplygin |63] as a convenient 
soluble model to study the lifting force on a plane wing 
in aerodynamics. This equation of state has also a con¬ 
nection with string theory and admits a supersymmetric 
generalization (see m for details). For the equation of 


9 This is most easily seen by considering Eq. |2j) with P = 0 and 
A > 0 leading to e m = 2eyv, hence e = 3 ca- 


state (21), the equation of continuity (|T|) can be inte¬ 
grated into 




(?)■“ 


1/2 


( 22 ) 


where a* is a constant of integration. Eq. ( |22| can be 
viewed as a combination of dark matter and dark energy. 
For a —> 0, we have e oc a -3 which behaves as pres¬ 
sureless dark matter. For a —► + 00 , we have e —> y[A 
which behaves as dark energy with a constant energy 
density e A = \[A , e quival ent to a cosmological constant 
(see Secs. IIC and IID). Although this model has an 


important historical interest (and can be motivated by 
string theory), it does not give a good agreement with 
the observations ED Eg. Therefore, generalizations of 
the Chaplygin gas model have been considered. 


F. Constant pressure 


As a first step, we consider a single dark fluid described 
by an equation of state corresponding to a constant neg¬ 
ative pressure 


P = -e A . (23) 

The equation of continuity ([!]) can be integrated into 


e = e A 


( t )- 1 


(24) 


where a* is a constant of integration. This model also 
provides a unification of dark matter and dark energy, 
and is simpler than the Chaplygin gas [34H37] . Actually, 
this model gives the same result as the ACDM model of 
Secs. IIC and II D| although its physical justification is 
different. As a result, this model agrees with the cos¬ 
mological observations, unlike the Chaplygin gas. In or¬ 
der to describe (small) deviations to the ACDM model, 
and in order to describe simultaneously the large scale 
structure of the universe and dark matter halos (as dis¬ 
cussed in the Introduction), we should consider equations 
of state that are close to, but different from, a constant 
pressure. This is the aim of the polytropic and logotropic 
equations of state considered in the following sections. 


G. Polytropic model of type I 

We assume that the universe is made of a single dark 
fluid. We can specify its equation of state in the form 
P = P(e), where e is the energy density. We consider a 
polytropic equation of state 

P = K &■ r = i + f W 

in which the pressure is a power of the energy density. 
This corresponds to the polytropic equation of state con¬ 
sidered in the first paper of Tooper [IT] in the context of 
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general relativistic stars. Therefore, we shall call it poly¬ 
tropic model of type I or, simply, Model I. This polytropic 
model has been studied in a cosmological context in US¬ 
ES- It can be viewed as a generalization of the Chaplygin 
gas model. 

For the equation of state ( [25] ), the equation of con¬ 
tinuity 0 can be integrated analytically. Assuming 
w = P/e > —1 (non-phantom), K < 0, and n < 0 (i.e. 
7 < 1 ), the relation between the energy density and the 
scale factor can be written as l35l: 


3/M 

— +1 


(t) 


i M 


(26) 


where e* = (|/\ \/c 2 )' n 'c 2 and a* is a constant of integra¬ 
tion. For n = —1/2, 7 = —1, and K = — A/c 2 


we recover 


the Chaplygin gas model (21) and (22). For n = — 1, 


7 = 0 , and I\ = —£a, we recover the constant pressure 
model (23) and (24), equivalent to the ACDM model. 


More generally, Eq. (26) can be viewed as a combination 


of dark matter and dark energy. Indeed, for a —> 0, we 
have e ~ e*(a*/a ) 3 which behaves as pressureless dark 
matter. For a —» + 00 , we have e —> e* which behaves 
as dark energy. This allows us to make the identification 
e* = £A- In this analogy, using Eq. (10) and the relation 


following Eq. (26), we find that the polytropic constant 


K is equivalent to a cosmological constant 


A = 8 ttG ( 


M 


(27) 


Introducing present-day variables, we can rewrite Eq. 


(26) as 


e = e A 


For a —> 0, 


, i/M 

^1 -1 

£A 



(28) 


e — £A 


e A ) 


i/M 


For a —¥ +00, 


- 1 


£A- 


M 


(?)’■ 


(29) 


(30) 


This model behaves as pressureless dark matter at early 
times and as dark energy at late times. However, this de¬ 
composition is only asymptotic and, in the intermediate 
regime, it is not possible to separate these two compo¬ 
nents unambiguously because they do not appear addi- 
tively in the energy density unless n = — 1 (i. e. 7 = 0 ), 
corresponding to the ACDM model (see Sec. IIF). We 
also note that the contribution of baryonic matter, con¬ 
sidered as a pressureless fluid (Pf, = 0), must be added 
in the total energy density as an independent species. 

The polytropic model (25) with K < 0 and n < 0 can 


be used to discuss deviations from the ACDM model (n = 


— 1). It is found from observations that n = —l-OS/lo ig 
l51l . If we impose that the velocity of sound c s is real, 
the condition c 2 = P'(e)c 2 = ( K'y/c 2 )(e/c 2 )' 1 ~ 1 > 0 con¬ 
strains the index n in the range — 1 < n < 0 (i.e. 7 < 0 ). 
Combining this condition with the observations, we get 
-1 < n < -0.9 (i.e. -0.11 < 7 < 0). 

In conclusion, the polytropic model of type I offers a 
unified picture of dark matter and dark energy [5ol 02] • 
In this model, the universe starts from the matter dom¬ 
inated epoch and approaches a de Sitter phase at late 
times. Contrary to the ACDM model, the cosmic coin¬ 
cidence problem (namely why the ratio of dark energy 
and dark matter densities is of order unity today) [Bl] is 
solved naturally in the polytropic gas scenario in E]. 
We finally note that a theoretical justification of a poly¬ 
tropic equation of state of the form of Eq. (25) for a per¬ 


fect relativistic fluid in cosmology has been given recently 
by Kazinski l65l in relation to the quantum gravitational 
anomaly. He predicts a polytropic equation of state of 
the form P = K{e/c 2 ) 1+1 ^ n with n € N. 


H. Polytropic model of type II 

As in the previous section, we assume that the uni¬ 
verse is made of a single dark fluid but, here, we specify 
its equation of state in the form P = P(p), where p is 
the rest-mass density. In order to solve the Friedmann 
equation 0 , we need to determine the energy density 
e. To that purpose, we assume that the dark fluid is at 
T = 0, or that the evolution is adiabatic (which is the 
case for a perfect fluid ED)- In that case, the first law of 
thermodynamics reduces to (see Appendix |A|): 




(31) 


For a given P(p) this equation can be solved to give the 
relation between the energy density and the rest-mass 
density. This relation can be written as (see Appendix 




e = pc 2 + p[ 1 - S V- dp' = pc 2 +u{p), (32) 


P' 


where pc 2 is the rest mass energy of the dark fluid and 
u(p) is its internal energy. On the other hand, combining 
the first law of thermodynamics (31) with the continuity 
equation ([!]), we get 


dp , „a 
— + 3 -p — 0 . 
at a 


(33) 


We note that this equation is exact for a fluid at T = 0, 
or for a perfect fluid, and that it does not depend on the 
explicit form of the equation of state P(p). It expresses 
the conservation of the rest-mass density p (or number 
density n = p/m). It can be integrated into 


( a oV 

P = P °{-J 


(34) 




























where po is the present value of the rest-mass density and 
ao is the present value of the scale factor. Substituting 


Eq. (34) in Eq. (32), we obtain 


. ( a o\ 3 | I" ( a o\ 

« = «c {-) +u «(-) 


(35) 


For any barotropic equation of state P(p), our treat¬ 
ment shows that the energy density in Eqs. (32) and (35) 
is the sum of two terms. The first term corresponds to 
the rest-mass energy of the fluid and the second term cor¬ 
responds to its internal energy. In the present approach, 
there is just one (relativistic) dark fluid. If we want to re¬ 
late our approach to the traditional viewpoint where the 
universe is made of dark matter and dark energy, we note 
that the first term in Eqs. ( f32j ) and (35) is equivalent to 
the term of Eq. © that describes a pressureless fluid 
(dark matter). As a result, we can interpret the second 
term in Eqs. ( [32] ) and ( |35[ ) as the term corresponding to 
the dark energy. The equation 


, 5 r p(p') ,, 

U (P) = P / 72 d P 

J p' 


(36) 


clearly shows that the dark energy term arises from pres¬ 
sure effects (i.e. collisions between particles). If P = 0, 
the internal energy vanishes and we recover the standard 
EdS mod el describing only pressureless dark matter (see 
Sec. IIBI. Alternatively, if we assume that the pressure 


is constant, P = — ca, we find that the internal energy is 
constant, u = e\, and we recover the model of Sec. |IIF| 
that is equivalent to the ACDM model. Therefore, the 
rest-mass energy pc 2 of the dark fluid mimics “dark mat¬ 
ter” and its internal energy u{p) mimics “dark energy” 
(see Appendix B of mi). This decomposition leads to 
a natural, and physical, unification of dark matter and 
dark energy, and elucidates their mysterious nature. We 
note that, in general, the dark energy term u{p) in Eq. 


(35) is not constant (it depends on the rest-mass den¬ 
sity p, and it does not even necessarily tend to a con¬ 
stant for p —> 0 ), c ontrary to the standard dark energy 


term e\ in Eq. (24) corresponding to a constant negative 
pressure. Therefore, our approach can be used to obtain 
generalized models of dark energy by considering differ¬ 
ent expressions of the equation of state P(p). If P(p) 
is sufficiently close to a constant, the resulting models 
will be relatively close to the ACDM model. Therefore, 
we can use this approach to study small deviations from 
the ACDM model and constrain the possible equations 
of state P(p) of dark energy from observations. 

To be specific, we consider a polytropic equation of 
state 


P = Kp\ 


1 

7 = 1 - 1 — 
n 


(37) 


in which the pressure is a power of the rest-mass density. 
This corresponds to the polytropic equation of state con¬ 
sidered in the second paper of Tooper [42] in the context 
of general relativistic stars. Therefore, we shall call it 


polytropic model of type II, or simply, Model II. This 
polytropic model has been studied in a cosmological con¬ 
text in [43} SI] , and independently in |43]. 


For the equation of state (37), using Eq. (32), we find 


that the energy density is related to the rest-mass density 
by 


= pc 2 + Kp In > (7 = 1), 


(38) 


e = pc 2 + — p y = pc 2 + nP(p), (7 ± 1). (39) 

7 - 1 


Combining Eq. (34) with Eqs. (38 1 and (39), we obtain 
for 7 = 1 : 


e = 


Poe 2 (-) 3 + (-) 3 In [- (-) 

\ a y V a / yp* \ a J 

and for 7 ^ 1 : 


e = Poe" ( — 


'(?)' 


K 7 / ao \ 37 


fPo . 

7 — 1 \ a 


er 


(40) 


(41) 


Introducing relevant (see below) notations Q m ,0 = 

p 0 c 2 /e 0 , n ' m)0 = Kp 0 /e 0 and ft 7 , 0 = ATpo /[(7 - !) e o], 
the Friedmann equation (|3| with A = 0 can be written 


as 


H_ _ A I »m ,0 | ^m,o 1 n[(Po/p*)(aoA0 3 ] 


H 0 \J (a/ao) 3 
for 7 = 1 , and as 

f- 


(a/a 0 y 


n 


m, 0 


n 


7,0 


(a/ao ) 3 (a/a 0 ) 37 


(43) 


for 7 ^ 1 . By construction, fi mi0 + 0 ln(po/p*) = 1 
and + 1 ^ 7,0 = 1. As discussed previously, the rest- 
mass energy of the dark fluid mimics dark matter so that 
fl m ,o represents the present fraction of dark matter. The 
contribution of baryonic matter, considered as a pres¬ 
sureless fluid (Pb = 0), must be added in the total energy 
density, as an independent species. Since it simply adds 
to the rest-mass energy of the dark fluid, we can take it 
into account by considering that n m ,o actually represents 
the present total fraction of matter (baryonic and dark). 

For the polytropic equation of state (37), the energy 
density in Eqs. (38)-(41) is the sum of two terms. An 
ordinary term e m = pc 1 oc a~ 3 (rest-mass energy) equiv¬ 
alent to pressureless dark matter and a new term e 7 = 
u(p) = P(p)/( 7 — 1) oc a -37 (internal energy) depending 
on the polytropic index 7. The equation of state of the 
new term is (see also Appendix [ b]) P = (7 — l)u = u/n. 
When 7 > 1 (i.e. n > 0), we find that P ~ e/n and 
e ~ nKp 1 oc a -37 in the “early” universe (small a, large 
p) and that P ~ A'(e/c 2 ) 7 and e ~ pc 2 oc a -3 in the 
“late” universe (large a, small p). When 7 < 1 (i.e. 
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n < 0), we find that P ~ A'(e/c 2 ) 7 and e ~ pc 2 * oc a~ 3 in 
the “early” universe (small a, large p ) and that P ~ e/n 
and e ~ nKp 7 oc a -37 in the “late” universe (large a, 
small p ). 10 Therefore, the new term (i.e. the internal 
energy of the dark fluid) dominates in the late universe 
when 7 < 1 (i.e. n < 0 ), while it dominates in the early 
universe when 7 > 1 (i.e. n > 0). This remark allows us 
to refine the discussion given previously. The new term 
can mimic dark energy only for 7 < 1 (i.e. n < 0). For 
7 = 0 (i.e. n = —1), the new term corresponds to a 
constant dark energy density (P = —u, u = ca), equiva¬ 
lent to the ACDM model. By contrast, when 7 > 1 (i.e. 
n > 0 ), the new term can mimic the effect of an exotic 
constituent appearing in the early universe. For 7 = 2 
(i.e. n = 1 ), the new term is equivalent to stiff matter 
(P = u, u = e s oc * a~ 6 ) as discussed in [33] in relation 
to self-interacting BECs at T = 0 and in relation to the 
cosmological model of Zel’dovich ( 66 , 67]. For 7 = 4/3 
(i.e. n = 3), the new term is equivalent to the radiation 
of an ultra-relativistic gas (P = u/ 3, u = e ra d oc a~ 4 ). 


We now restrict ourselves to the case 7 < 1 in order to 
describe dark energy. If we assume that the energy den¬ 
sity of the new term (internal energy) is positive in order 
to account for the observations (some energy is missing 
with respect of the EdS model where there is only dark 
matter), we conclude from the relation e 7 = P/(y — 1) 
that the pressure must be negative (K < 0, P < 0). 
This is the case, in particular, for the constant equation 
of state P = —£a < 0, corresponding to 7 = 0, equiva¬ 
lent to the usual ACDM model. If we assume that the 
dark fluid is made of a self-interacting scalar field (e.g. a 
BEC), it is shown in Appendix [C| that a negative pres¬ 
sure can arise from an attractive potential of interaction. 
Therefore, an equation of state with a negative pressure 
(such as the equation of state of dark energy) can be 
justified from a field theory. 


The polytropic model (411 with an index 7 < 1 can be 
used to discuss deviations from the ACDM model (7 = 
0). A detailed study has been done recently in [35] ■ It is 
found from cosmological constraints that —0.089 < 7 < 0 
(corresponding to —1 < n < —0.92). 11 We note that this 
range is very close to the one obtained from Model I (see 
Sec. IIG) but this may be a coincidence. 


Before closing this discussion, we briefly recall the 
analogies and the differences between the two polytropic 
models: (i) Models I and II describe dark energy provided 
that 7 < 1 (i.e. n < 0) and A < 0 (negative pressure). 
They are equivalent to the ACDM model when 7 = 0 
(i.e. n = —1). (ii) In Model II, the energy density (41) 


10 Since P ~ (7 — 1 )e in the late universe, we find that w = P/e —> 

7 — 1 when a —> +oo. Therefore, the universe is normal when 

7 > 0 (the energy density decreases as the scale factor increases) 

and phantom when 7 < 0 (the energy density increases as the 

scale factor increases). 

11 For 7 < 0, the universe becomes phantom at late times (see 

footnote 10 ). 


is the sum of two terms (rest-mass energy and internal 
energy) that can be associated with what has been called 
dark matter and dark energy, respectively. In Model I, 
this distinction is less clear because the expression (28) 
of the energy density is non-additive, so the distinction 
between dark matter and dark energy is only asymptotic, 
(iii) In model I, the energy density tends to a constant ca 
at late times (leading to a de Sitter era), which is equiva¬ 
lent to the cosmological constant or to the ordinary form 
of dark energy. In model II, the dark energy term is not 
constant (i.e. it depends on time), even asymptotically. 


I. Cardassian model 


There is an interesting consequence of the ideas devel¬ 
oped in the previous section. Using Eq. (32), we can 
write the Friedmann equation ([3]) with A = 0 in terms of 
the rest-mass density p as 


tt2 87 tG 

H = —z~p + v(p). 


(44) 


where p = po(ao/a) 3 and v{p) = ( 8nG/3c 2 )u(p ). This 
equation is similar to the modified Friedmann equation 
introduced by Freese and Lewis [58] in the so-called Car¬ 
dassian model . 12 For the polytropic equation of state 
(37), it exactly concides with the power law Cardassian 


model. Eq. (44) was originally [ 68 ] justified by the fact 
that the Friedmann equation is modified as a consequence 
of embedding our universe as a three-dimensional sur¬ 
face (3-brane) in higher dimensions. Hence, the modified 
Friedmann equation (44) may result from the existence of 


extra-dimensions. Our approach provides another, sim¬ 
pler, justification of this equation from the ordinary four 
dimensional Einstein equations. In particular, starting 
from the usual Friedmann equation <© and considering 
a dark fluid at T = 0, or an adiabatic fluid, the “new” 


term v{p) in the “modified” Friedmann equation (44) can 
be interpreted as the internal energy (36) of the dark fluid 


while the “ordinary” term ( 8irG/3)p corresponds to its 
rest-mass density. 


III. LOGOTROPIC EQUATION OF STATE 


In Sec. IIH we have argued that dark matter and dark 
energy may be the manifestation of a unique dark fluid 
characterized by an equation of state P(p). If this claim 
is correct, the equation of state P(p) should describe both 
the cosmological evolution of the universe and the struc¬ 
ture of the dark matter halos. Let us therefore consider 
the possibility that dark matter halos are described by an 
equation of state P(p), the same as the one governing the 


12 The name Cardassian refers to a humanoid race in Star Trek 
whose goal is to take over the universe m ■ 
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cosmological evolution of the universe. In that case, they 
satisfy the classical condition of hydrostatic equilibrium 
between pressure and gravity 13 

VP + pV<F = 0 . (45) 


We have seen that a dark fluid with a constant pressure 
(P = — e\) generates a cosmological model equivalent 
to the ACDM model. Therefore, this equation of state 
seems to be suitable at the cosmological scale. However, 
when this equation of state is applied to dark matter 
halos, we find that VP = 0, so there is no pressure gra¬ 
dient to balance the gravitational force. Therefore, as in 
the classical CDM model where P = 0, the dark matter 
halos described by a constant equation of state display 
cuspy density profiles while observations tend to favor a 
constant density core (cusp/core problem). This obser¬ 
vational result is a strong argument against a constant 
equation of state, thus against the ACDM model, be¬ 
cause it is not suitable at the scale of dark matter halos. 
Therefore, we turn next to a polytropic equation of state 
of the form of Eq. (371. The condition of hydrostatic 
equilibrium (45) becomes 


temperature. We note that A must be positive if we 
want that the pressure forces balance the gravitational 
attraction. This is satisfying from a thermodynamical 
point of view. In the following sections, we successively 
apply the logotropic equation of state to the evolution of 
the universe and to the structure of dark matter halos. 


IV. LOGOTROPIC COSMOLOGY 

A. The equation of state 

We assume that the universe is made of a single dark 
fluid described by a logotropic equation of state 


P = A In 


(49) 


where p is the rest-mass density, A is the logotropic tem¬ 
perature (we assume A > 0, see Sec. Ill I, and p* is a 
reference density . 15 It will be called the logotropic dark 
fluid (LDF). According to Eq. (32), the relation between 
the energy density and the rest-mass density is 


ATyp 7 1 Vp + pV'I> = 0 . 


(46) 


Because of pressure gradients when If / 0 and 7 ^ 0, 
the dark matter halos described by a polytropic equation 
of state display a central core instead of a cusp . 14 If we 
want this equation of state to describe both dark matter 
halos and the cosmological evolution of the universe, 7 
must be close to zero (i.e. the equation of state must 
be close to a constant pressure, equivalent to the ACDM 
model). But, in that case, we encounter the problem of 
the justification. Indeed, it seems difficult to theoreti¬ 
cally justify why 7 should take a value like —0.089 for 
example [45]. Therefore, we consider another possibility. 
We assume that 7 —> 0 and K —> 00 in such a way that 
A = Kj is finite [53]. In that limit, Eq. (46) becomes 


Vp 

A— + pV<F = 0 . 
P 


(47) 


Comparing Eq. (47) with Eq. (45), we see that the 
pressure term corresponds to the logotropic equation of 
state 


P = Alnp + C, (48) 

where A and C are constants. As shown in Appendix 


C 5 c the parameter A can be interpreted as a logotropic 


13 For dark matter halos, we can use Newtonian gravity. 

14 For the pressure gradient to counterbalance the gravitational at¬ 
traction, we need if 7 > 0. Since we have shown from cosmolog¬ 
ical considerations in Sec. |II H| that K < 0, we must have 7 < 0. 
This is consistent with the cosmological bound 7 < 0 obtained 
in 1451 . 


The pressure is related to the internal energy by P = 
—u — A. Using Eq. (341, we obtain 


£ = Poc 2 —Ain 


'po 

( a o \ 3 

p* 

\ a ) 


- A. (51) 


As explained in Sec. |II H[ the first term (rest-mass en¬ 
ergy) in Eqs. (501 and (|51[) mimics dark matter and the 


second term (internal energy) mimics dark energy. The 
evolution of the energy density with the scale factor is 
plotted in Fig. [l]and the relation between the energy den¬ 
sity and the rest-mass density is plotted in Fig. [2] The 
universe starts at a = 0 with an infinite rest-mass den¬ 
sity (p —» +00) and an infinite energy density (e —> +00). 
The rest-mass density decreases as a increases, see Eq. 
( [34| ). The energy density decreases as a increases (i.e. p 
decreases), reaches a minimum cm = — A ln(A/p*c 2 ) at 
«m = a 0 (poc 2 /A) 1//3 (i.e. pu = A/c 2 ), increases as a in¬ 
creases (i.e. p decreases) further, and tends to e -> +00 
as a — > +00 (i.e. p — > 0). The branch a < «m (i.e. 
P > Pm) corresponds to a normal behavior in which the 
energy density decreases as the scale factor increases. 
The branch a > cim (be. p < Pm) corresponds to a 
phantom behavior in which the energy density increases 
as the scale factor increases. We note that A is equal 
to the rest-mass energy at the point where the universe 
becomes phantom. 


15 The meaning of this density will be unraveled in Sec. |Vl| but we 
leave it unspecified for the moment. 
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Combining Eqs. (49) and (501, we obtain the following 


relation between the pressure and the energy density 


e = p 


,c 2 e 


P/A - P - A. 


(52) 


This determines the equation of state P = P(e) under 
the form e = e(P). The pressure decreases as a increases 
(i.e. p decreases). It starts from P —► Too at a = 0 (i.e. 
p — > +oo, e — > +oo), achieves the value Pm = -£m at a im 
(i.e. pm, £m) and tends to P —> —oo when a —► +oo (i.e. 
p —> 0, e —> Too), see Fig. [4j The equation of state P(e) 
is defined for e > cm and has two branches corresponding 
to a normal universe (P > Pm) and a phantom universe 
(P< Pm), as shown in Fig. [5] Therefore, the equation 
of state P(e) is multi-valued. 

In the early universe (a —» 0, p —> Too), the rest-mass 
energy (dark matter) dominates, and we have 


2 2 / a o V 

e~pc ~ poc — j 


P ~ Ain 


P*c 


(53) 


For very small values of the scale factor, we recover the 
results of the CDM model (P = 0) since e oc a -3 . In 
the late universe (a -T Too, p —> 0), the internal energy 
(dark energy) dominates, and we have 


-Ain 


P* 


3 A In a, 


P~ -e. 


(54) 


We note that the equation of state P(e) reduces to P~ 
—e for e —> Too, which is similar to the usual equation 
of state (111 of the dark energy . 16 However, this is only 


an asymptotic result, different from the usual equation of 
state of the dark energy where P = — e exactly, implying 
a constant energy density e = e A (see Sec. IID). As a 


result, in our model, the energy density does not tend to 
a constant e A for a —> Too, but slowly increases as In a. 
This corresponds to a phantom regime where P/e < — 1 
at high energy densities (we have PTe—t—A<0 when 
e -T Too on the phantom branch). As we shall see, this 
leads to a super de Sitter behavior. We also emphasize 
that a polytropic equation of state with an index 7 = 0 
(or n = —1), corresponding to a constant pressure P = 
—e A , leads to different results since the energy d ensit y 
tends to a constant value e A for a —> Too (see Sec. IIF I, 


which is different from Eq. (541. 


B. The energy density 

Since, in o ur m odel, the rest-mass energy of the dark 
fluid in Eq. (511 mimics dark matter, we identify p 0 c 2 


with the present energy density of dark matter. We set 
Poc 2 = flm.oeo, (55) 


where eo is the present energy density of the universe and 
n mfi is the present fraction of dark matter . 17 As a re¬ 
sult, the present internal energy of the dark fluid uq = 
£o — Poc 2 is identified with the present dark energy density 
e A = (1 — Sl mj o)eo where fl Aj o = 1 — fi m ,o is the present 
fraction of dark energy. From the observations, one has 


H 0 = 70.2kms _1 Mpc -1 = 2.27510“ 18 s" 1 , 0 = 

r2„2, 


0.274 and f2 Aj0 = 0.726, yielding eo = 3H§c 2 /8irG = 
8.32 x 10~ 7 gm - 1 s -2 and cq/c 2 = 9.26 x 10 - 24 gm -3 . 
Therefore, po = 2.54 x 10 - 24 gm -3 and p A = e A /c 2 = 
6.72 x 10 - 24 gm“ 3 . 


Applying Eq. (51) with Eq. (551 at a = a 0 , we get 


Ti — ^m,o£o — ^4 In 


^m.CTO 
P*C 2 


-A, 


(56) 


which determines /?* as a function of A, for known (mea¬ 
sured) values of eo and £l m ,o- Using Eqs. (55) and (56), 


the relation (|5l|) between the energy density and the scale 

3A 


factor can be rewritten as 
c _ U m o 


{a/a o y 


T (1 — U mj o) 


IT 


eo(l — U m0 ) \ a o 


In I — 


(57) 

It is convenient to introduce the dimensionless logotropic 
temperature 


P=A = 


£o(l U m o) 
and the normalized scale factor 


R=—. 

a 0 


(58) 


(59) 


In terms of these variables, the characteristic density p* 
is given from Eq. (56) by 


p*c 




= e 1+1 / B , 


m, 0 


(60) 


and Eq. (57) becomes 
£ _ U m o 


£0 


P 3 


T (1 — U m , 0 )(l T 3P In R). (61) 


The dimensionless logotropic temperature B is the only 
free parameter of the model. 

For B = 0, Eq. (61) reduces to 


a 


m, 0 


eo 


R 3 


T 1 - fi 


m, 0 • 


(62) 


This returns the ACDM model with a different justifica¬ 
tion. This is because when B —>• 0, using Eqs. ((581) and 


16 It is interesting to note that the equation of state P = — e can 
be ob tain ed as an asymptotic limit of the logotropic equation of 
state |49|). This was not apparent at first sights. 


17 As explained in Sec. |ll H[ we also include baryonic matter in 
po so that actually represents the present total fraction of 

mass (baryonic and dark). 
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FIG. 1: Evolution of the energy density as a function of the 
scale factor. 



pc7e 0 

FIG. 2: Relation between the energy density and the rest- 
mass density (for illustration we have taken B = 0.2). 


(60), the equation of state (49) reduces to a constant 
negative pressure P = — A/B = —eo(l — f2 m ,o) = — £a 


equivalent to Eq. (23). 
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For B ^ 0, our model differs from the ACDM model. 
The relation between the energy density and the scale 
factor is represented in Fig. [I] for different values of B. 
For R —> 0, 


e fl 


m,0 


eo R 3 


(64) 


and, for R —>• +oo, 


-~3B(l-fi rni0 )lni2, (B^O). (65) 

eo 

When B > 0, the curve e(R) presents a minimum at 
{Rm,£m), as detailed in the next section. When R < 
Rm, the energy density decreases with the scale factor, 
which corresponds to a “normal” universe. When R > 
Rm, the energy density increases with the scale factor, 
which corresponds to a “phantom” universe. 

It is useful to write the relation between the energy 
density and the rest-mass density using dimensionless 
variables. According to Eqs. (34), (55) and (59), we 
have 


0^0 I 

P = -2 

c z 


R 3 ' 


( 66 ) 


Eqs. (61) and (66) determine the relation e(p) in para¬ 
metric form. Eliminating R between these two relations, 


18 More precisely, using Eqs. (|58| and ( |60| ), we can write the lo- 
gotropic equation of state (|49|) under the form 


P = Be a In p - e A 


1 + B + B In 


f eo H77 


(63) 


For B = 0, we get P = —ea which is equivalent to the ACDM 
model (see Sec. II Fi. 


we explicitly obtain 


eo 


eo 


+ (1 — 


1 + B In 


^m.0^0 


pc* 


(67) 


as represented in Fig. [2] For B = 0, this equation reduces 
to e = pc 2 + ca equivalent to the ACDM model. 


C. The point of minimum energy density 


The minimum of the curve e(i?) defined by Eq. (61) is 
given by 

1/3 


Rm — 


— ) — (1 —f!m,o) 

e 0 / M 


n 


m, 0 


_B{ 1 — 

B + 1 + B In 


a 


m, 0 


1 - n 


m,0 


( 68 ) 


-BlnB 


( 69 ) 

The function €m(B) is represented in Fig. |3] For the 
ACDM model, corresponding to B = 0, we have 


(70) 


- (0) = 1 - n m , o = 0.726. 

£ 0 J M 


In that case, Rm +oo, so the minimum cm = Gv is 
rejected at infinity. For B —> 0, 


- = (1 — ^m,o)(l — BlnB), (71) 

£ 0 J M 


and for B —» +oo, 


eo 


= — (1 — Q m> o)-Bln.B — > — oo. (72) 


M 


We note that the function cm{B) first increases with B 
before decreasing. Its maximum is located at 


B, = 


n 


771,0 


i - a 


771,0 


= 0.377, 


- (B0 = 1. (73) 

e 0 J M 
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The function cm{B) recovers its “initial” value cm( 0) at 
the point 


B 2 = 




m, 0 


i - n 


-e = 1.03, 


m, 0 


(74) 


- (P 2 ) = 1 - = 0.726. (75) 

<+/ M 


Finally, the minimum energy density cm(P) vanishes at 
Bmajf given by the implicit equation 

ln(-Bmax) - * = 1 + In f ^ . (76) 


B „ 


1 - a 


m,0 


Solving this equation numerically, we find P max = 1.79. 
In the following, we assume 0 < B < B max = 1-79 so the 
energy density e(R) remains always positive (which is of 
course a necessary condition). This puts a first constraint 
on the allowable values of B. If we demand that the 
present-day universe is not phantom, corresponding to 
the condition R M > 1, using Eqs. (l68|) and ([73|, we 


obtain the more stringent constraint 0 < B < B\ = 
0.377. This gives a physical interpretation to B\. 



FIG. 3: Minimum energy density tM (B) as a function of the 
logotropic temperature B. 


D. The pressure 


Using Eqs. (34), (55), (58), (59) and (60), the pressure 


(49) can be expressed in terms of the scale factor as 

P = —e 0 (l - n mi0 )(B + 1 + 3 B Ini?). (77) 


The pressure vanishes at 

Rw = e -( B+1)/3B . (79) 


At that point, according to Eq. (61), the energy density 
is 


(-) = u m , oe ( B+1 >/ B - (1 -fi ro , 0 )B. (80) 

The pressure is positive for R < R w and negative for 
R > R w . Since R w < 1, the pressure is always negative 
in the present-day universe, for any B > 0. Its value is 
Po = — (B + 1)ca- The relation between the pressure and 
the scale factor is plotted in Fig. [4] for different values of 
B. We note the “exceptional” point R e = e -1 / 3 = 0.7165 
at which the pressure takes the same value P e = — e a = 
—6.04 x 10 -7 gm -1 s -2 , independently of B. 



R 


FIG. 4: Evolution of the pressure as a function of the scale 
factor. 



For B = 0, we obtain a constant pressure P = — e 0 (l — 
fl m , 0 ) = — £a, equivalent to the ACDM model. For B ^ 
0, our model differs from the ACDM model. For R —> 0 
and for R -+ +oo, 

(78) 


FIG. 5: Equation of state P(e) giving the pressure as a func¬ 
tion of the energy density (for illustration we have taken 
B = 0.2). 


The equation of state P(e) is given in parametric form 
ry Eqs. (61) and ( [77| . Using Eqs. (52), (58) and (60), 


P 


3Pe 0 (l - f2 m , 0 ) In R. 
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we also have 


— = Q, m oe (B + l)/B e P/[e 0 (l-n m , 0 )B] 

e 0 


(1 ^m,o) 


P 

£o(l ^m,o) 


+ B 


(81) 


which determines e(P) in dimensionless form. The equa¬ 
tion of state P(e) presents two branches corresponding 
to a normal universe and a phantom universe (see Fig. 
[5]). For e —¥ +00, we get 

P ~ Pe 0 (l - fi m , 0 ) lne (82) 

on the normal branch (P —> 0) and 

P ~ -e (83) 

on the phantom branch (P —► +00). The pressure at the 
point of minimum energy 6m is 


Pm — £o(l 0) 


B + 1 + B In 


n, 


m,0 


l-n 


m, 0 


-PlnP 
(84) 

Comparing E qs. fl69| ) and ([84]), we find that Pm = — cm 
( see also Sec. IVA|. 


E. The equation of state parameter w 


The equation of state parameter w is defined by 


P = we. 


(85) 


According to Eqs. (61) and {77} , it can be expressed in 
terms of the scale factor as 


— (1 — fl m 0 )(P + 1 + 3PlnP) 

w = ~o -• (86) 

%^ + (l-fi ro , 0 )(l + 3PlnP) 

For P —> 0, 

w ~ -3P 1 ~° m ’° p 3 hiP. (87) 

^ 'm,0 

For P —>• +00, 

w ->• -1. (88) 

More precisely, w + 1 ~ —1/(3 In P) for P —> +00. We 
note that w > 0 for P < R w and w < 0 for P > R w . 
Therefore, the pressure is positive for P < R w and neg¬ 
ative for P > R w , in agreement with the results of Sec. 
|IVD[ We also note that w > — 1 for P < Rm and 
w < — 1 for P > Rm- Therefore, the universe is nor¬ 
mal for P < R m (the energy density decreases as the 
scale factor increases) and phantom for P > Rm (the en¬ 
ergy increases as the scale factor increases), in agreement 
with the results of Sec. IIVBI 

The curve w(R) is plotted in Fig. [6j It starts from 
w = 0 at P = 0, increases, reaches a maximum (see Fig. 



FIG. 6: The curve w(R) for B = 0, B = 0.1, and B = 0.2 
(top to bottom). 



FIG. 7: Zoom of Fig. [6]close to R = 0 showing the maximum 
of the function w(R). 


[7|, decreases, becomes negative after P = R w , passes 
below —1 after P = Rm-i reaches a minimum, increases, 
and tends to —1 as P —>• +00. 

For B = 0, corresponding to the ACDM model, the 
function w(R) is given by 


(1 Hm,o) 

+1 - n 


R 3 


m, 0 


(89) 


It starts from w = 0 at P = 0, decreases monotonically, 
and tends to — 1 as P —> +00. 

The present value of the equation of state parameter 
is 


Wo = -(l-fi mi0 )(P + l). (90) 

For B = 0, corresponding to the ACDM model, Wq = 
-(l-fim,o) = -0.726. 
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F. The velocity of sound 


The velocity of sound c s is defined by c 2 = P'(e)c 2 . 
Taking the derivative of Eq. (52) with respect to e, we 
obtain 


P£l _ 1 
A 1 


(91) 


This relation requires that A > 0, hence B > 0, otherwise 
the velocity of sound would always be imaginary (c 2 < 0). 
The constraint A > 0 is consistent with the results of 
Sec. |III| It also reinforces the interpretation of A as a 
(logotropic) temperature. 19 Substituting Eqs. (34), (58) 
and (|68| in Eq. (91), we get 


1 


c 2 (*#.)*-i 


(92) 


For B = 0, corresponding to the ACDM model, c 2 = 0 
since the pressure is constant (P = —e\). For B > 0, 
the velocity of sound is real for R < Rm (he. when 
the universe is normal) and imaginary for R > Rm (he. 
when the universe is phantom). The relation between the 
velocity of sound and the scale factor is plotted in Fig. 

m 



If we require that (c 2 )o > 0, we must have Rm > 1, 
hence B < Bi = 0.377. The condition that the present- 
day velocity of sound is real coincides with the condition 
that the present-day universe is non-phantom. 

More stringent constraints on B can be obtained from 
the following arguments. In order to have (c 2 /c 2 )o < 1, 
the dimensionless logotropic temperature must satisfy 


B < 


n 


m, 0 


2(1 — n mi0 ) 


Pi 

2 


0.1885. 


(94) 


This ensures that the velocity of sound in the present 
universe is less than the speed of light. In order to have 
(c 2 /c 2 )o < 1/3, the dimensionless logotropic temperature 
must satisfy 


B < 


n 


m, 0 


4(i - n m , 0 ) 


— = 0.09425. 
4 


(95) 


This ensures that the present universe is non-relativistic 
(i.e. the velocity of sound in the present universe is less 
than the velocity of sound in a radiation-dominated uni¬ 
verse for which P = e/3 and c 2 /c 2 = 1/3). We shall 
take the value Pi/4 = 0.09425 as an upper bound for 
the allowable values of B. 

For a given P, the velocity of sound becomes larger 
than the speed of light (c s > c) when R > R s with 


Rs 


a 


m, 0 


1/3 


2P(1 - On, 0 ) 


Rm 
2V3' 


(96) 


Since R s < Rm, the logotropic model may break down 
before the universe reaches the phantom regime at R = 
Rm- 

Remark: If we ompute the velocity of sound from the 
relation c 2 = P'(p), we obtain 


c 


2 

s 


A 

P 



(97) 


This approximation is valid when p A or R <C Rm ■ 


G. The deceleration parameter 


FIG. 8: Evolution of the square of the velocity of sound as 
a function of the scale factor (for illustration we have taken 
B = 0.2). 


The velocity of sound in the present-day universe is 

1 

E>3 

0 


K- 1 ' 


(93) 


In a flat universe without cosmological constant (k = 
A = 0), the deceleration parameter q = —da/a 2 is given 
by (see Eqs. |2]), ([3| and (85)): 


9 = 


1 + 3u> 


(98) 


Using the expression of w given by Eq. (86), we obtain 


s ^-(l-f] m , 0 )(3P + 2 + 6PlnP) 

9=-nr^-—• (99) 


^ + (l-n m , 0 )(l + 3Plnfl) 


19 For a classical isothermal gas with P = phsT/m and c^ = 
P'(p) = ksT/m , the velocity of sound is real at positive temper¬ 
atures and imaginary at negative temperatures, implying that 
negative temperatures are (usually) forbidden. 


For R —> 0, 


1 - 9 B 1 n nm ’° R 3 In R 


n 


m, 0 


( 100 ) 
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For R —> +oo, 


<?->-!• 


( 101 ) 


The deceleration parameter q(R) vanishes at R = R c de¬ 
termined in the next section. When R < R c , the universe 
is decelerating (q > 0) and when R > R c , the universe is 
accelerating (q < 0). 

The curve q(R) is plotted in Fig. |9| Its b ehavior fol¬ 
lows the one of w(R) described in Sec. |IVE 

For B = 0, corresponding to the ACDM model, the 
function q(R) is given by 


q = 


R 3 


— 2(1 — n m>0 ) 


fl 3 


i - a 


m,0 


( 102 ) 



FIG. 9: Evolution of the deceleration parameter as a function 
of the scale factor for B = 0, B = 0.1, and B = 0.2 (top to 
bottom). 


the value of B. We also note that, in the framework 
of the logotropic model, R c lies in the small interval 
[0.574, 0.6065]. Since R C (B) increases, we note that when 
B > 0 the acceleration of the universe occurs later than 
in the ACDM model. 



FIG. 10: The scale factor R c at which the universe accelerates 
as a function of B. 


The present value of the deceleration parameter is 

^m,o — (1 — ^m,o)(3-B + 2) 
qo = - 2 -’ V 1 U 6 ) 

For B = 0, corresponding to the ACDM model, 

q 0 = 3Slm ’°~ 2 = -0.589. (107) 

The condition for the present-day universe to be accel¬ 
erating is that Q m , 0 < 2/3. This inequality is indeed 
realized by the observational value = 0.274. 


H. The point at which the universe accelerates 


According to Eq. (99), the point R C {B) at which the 
universe starts accelerating (q = 0) is determined by the 
equation 


B = 


1 


-2 


3(1 + 21ni? c ) ‘ 

For B = 0, corresponding to the ACDM model, 

1/3 


R c = 

For B —»• +oo, 


a 


m, 0 


R 


2(1 — fl m> o)_ 

e" 1 / 2 = 0.6065. 


= 0.574. 


(103) 


(104) 


(105) 


The curve R C {B) is plotted in Fig. 10 Since R c < 1, 


the present-day universe is always accelerating, whatever 


I. The transition between matter and dark energy 


If we interpret in Eq. (61) the rest-mass density as 


“dark matter” and the internal energy as “dark energy”, 
their ratio evolves with the scale factor as 


DE 

DM 


u(p) 

pc 2 


i - n 




^i? 3 (l + 3Slni?). (108) 


m, 0 


The transition between dark matter and dark energy 
takes place at a scale factor R? determined by 


a 


^ = (l-n m , 0 )(l + 3Bln/? 2 ). 


(109) 


For B = 0, corresponding to the ACDM model, we 
get R 2 (0) = 0.723. For B —> +00, we get R 2 ~ 
exp {(2f2„^ 0 — 1)/3(1 — n mj0 )-B} — > 1. More generally, 
the function R 2 {B) is plotted in Fig. 11 
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FIG. 11: The scale factor R 2 corresponding to the transition 
between dark matter and dark energy as a function of B. 



FIG. 12: Temporal evolution of the scale factor for B = 0, 
B = 0.1, and B = 0.2 (bottom to top). 


J. The temporal evolution of the scale factor 


K. The age of the universe 


Using Eq. (61), the Friedmann equation <© with A = 0 
takes the form 


Taking R = 1 in Eq. (Ill), the age of the universe is 
given as a function of the logotropic temperature by 


H = I = H 0 ^^f + {l-tt mfi )(l + W\nR). (110) 

The temporal evolution of the scale factor R(t) is given 
by 


dx 


^ + (i - U mi0 )(l + SBlnx) 


= H 0 t. (Ill) 


For t —» 0, 



e 

eo 


4 


( 112 ) 


like in the EdS universe (see Sec. 


IIB). 


For t —» + 00 , 


to(B) — — 


nRo = l 


dx 


+ (1 — r2 m? o)(l + 35 111 x ) 
(114) 



(H3) 

This solution, which is valid in the regime where the uni¬ 
verse is phantom, has a super-de Sitter behavior. How¬ 
ever, as indicated in Sec. |IVF[ the logotropic model may 
break down before the universe enters in this regime be¬ 
cause the velocity of sound exceeds the speed of light 
when R> R s with R s < Rm- 


‘d>B 

-^-(1 — 


The curve R(t) is plotted in Fig. 12 For B = 0, corre¬ 
sponding to the ACDM model, we recover the analytical 
solution of Eq. (fl4|) . 


FIG. 13: The characteristic times of the logotropic model as 
a function of B (see the text for details). 

For B = 0, corresponding to the ACDM model, we 
recover Eq. leading to the standard value t o (0) = 

13.8 Gyrs. More generally, the function to(B) is plotted 
in Fig. [l3j We also plot the functions t w (B), t c (B), 
t 2 (B), tJJ3) an d tAf(S) corresponding to the time at 
which (i) the pressure becomes negative, (ii) the universe 
accelerates, (iii) the transition from dark matter to dark 
energy occurs, (iv) the velocity of sound becomes higher 
than the speed of light, (v) the universe enters in the 
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phantom regime. They are obtained by replacing Rq = 
1 in Eq. (Ill) by R Wl R c , R 2 , R s and Rm given by 
Eqs. ([79]), ( |103| ), (fl09), (96|), and ([68|). For B = 0, 


we recover the standard values i c (0) = 7.18 Gyrs and 
< 2 ( 0 ) = 9.61 Gyrs of the ACDM model. For B —> 0, we 
have the asymptotic behavior 


t w (B) 




' 2B — } 0 . 


771,0 


(115) 


L. The Hubble diagram 

In terms of the redshift parameter 



1 — ^0 
a 


the Hubble function (110) can be written as 


(116) FIG. 15: Hubble (/x versus z) diagram corresponding to the 
logotropic model with B = 0, B = 0.1 and B = 0.2 (bottom 
to top) compared with observational data. 


H(z) = H 0 ^n m , 0 (z + l) 3 + (1 - fi m , 0 )[l - 3Bln(z + 1)]. 

(117) 

The Hubble function H(z) is plotted in Fig. 14 for dif¬ 
ferent values of B and compared with the observational 
data of Refs. [53 Eg. 



source is then given by [72]: 




5 log 


Mpc 


+ 25. 


(119) 


The function /x(z) corresponding to the logotropic model 
is plotted in Fig. [U)]for different values of B , and com¬ 
pared with observational data of the Union 2.1 Compila¬ 
tion m- Even for values of B as large as 0.2, the curves 
are very close to the ACDM model and can hardly be 
distinguished. They are all consistent with the observa¬ 
tional results. 


M. The shift parameter from CMB 

The shift parameter R, which is related to the position 
of the first acoustic peak in the power spectrum of the 
temperature anisotropies d, provides an information 
from the CMB [45]. It is defined by 


FIG. 14: The Hubble function H(z ) corresponding to the 
logotropic model with B = 0, B = 0.1 and B = 0.2 (top to 
bottom) compared with observational data. 


The history of the universe expansion is revealed by 
the Hubble diagram giving the distance modulus /x(z) = 
to — M as a function of the redshift z, where to. is the 
apparent luminosity and M the absolute luminosity of 
a light-emitting source (standard candle) such as super¬ 
novae of Type la. The procedure to obtain this curve is 
explained in Ref. [45]. In a spatially flat universe, the 
luminosity distance of a source with redshift z is [71]: 


Mpc 


c(l + z) [ 
Jo 


dz' 

W)' 


(118) 


The iv-corrected distance modulus of a light-emitting 


R — \/ ^m, 0 


dz 


H(z)/H o’ 


( 120 ) 


where z* is the value of the cosmological redshift at pho¬ 
ton decoupling. We adopt the nine-year WMAP survey 
final result z* = 1091.64+0.47 TSi. The value of the shift 
parameter obtained by these authors is i? 0 bs = 1.7329 ± 
0.0058. The ACDM model gives Racdm = 1-7342. On 
the other hand, using a combination of the Planck first- 
data release with those of the WMAP survey El, the 
value R' ohs = 1.7407 ± 0.0091 is obtained. The curve 
R(B) obtained from the logotropic model is plotted in 
Fig. [Tg] From the observational value R a b s we find the 
constraint B < 0.0262, and from the observational value 
R'o bs we find the constraint B < 0.0379. These values are 
about 3 times smaller than the bound Hi/4 = 0.09425 
obtained in Sec. IIVFI 
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FIG. 16: The CMB shift parameter R as a function of the 
logotropic temperature B. The horizontal lines correspond to 
the upper bound R 0 b a = 1.7387 of the R -value obtained from 
the nine-year WMAP survey final results, and to the value 
R'obs = 1-7407 arising from the combination of the Planck 
first-data release with those of the WMAP survey. 


where a is a constant. 20 On the other hand, Eq. 
takes the form 


(123) 


A 



4t tG 


A 


P- 


(125) 


Assuming spherical symmetry, and introducing the nota¬ 
tions 


6 = 


Po 

•> 

P 



where po is the central density and 


r o = 


A 

AnGpo 


1/2 


is a typical core radius, we obtain 


1 d / 2 d0\ 1 

V d€J = 0 ’ 


with 


(126) 


(127) 


(128) 


V. LOGOTROPIC DARK MATTER HALOS 


0(0) = 1 , 0'(O) = 0. (129) 


In the previous section, we have described the evo¬ 
lution of the universe as a whole, assuming that it is 
made of a single dark fluid with a logotropic equation 
of state. We have argued that this equation of state de¬ 
scribes simultaneously dark matter and dark energy. In 
this section, we consider dark matter halos which form 
during the nonlinear regime following the gravitational 
instability (Jeans instability) of the uniform cosmologi¬ 
cal background. We consider the possibility that dark 
matter halos are also described by the logotropic equa¬ 
tion of state (491, with the same value of the logotropic 


temperature A, and we investigate the consequences of 
this assumption. 

For dark matter halos, we can use Newtonian gravity. 
Combining the condition of hydrostatic equilibrium 


VP + pV$ = 0 
with the Poisson equation 

A'b = 4-7T Gp, 


( 121 ) 


( 122 ) 


we obtain the differential equation 


V • = -4nGp. 


(123) 


For the logotropic equation of state (49), using Eq. (121) 
we find that the density is related to the gravitationa 
potential by the Lorentzian-type distribution 


P( r) = 


a+ i$(r)’ 


(124) 


This equation coincides with the Lane-Emden equation 
of index n = — 1 m- We note that the Lane-Emden 
equation of index n = — 1 cannot be obtained from a 
polytropic equation of state of index n = — 1 (which cor¬ 
responds to a constant pressure P = K ) because, when P 
is a constant, there is no pressure gradient in Eq. |m| ) to 
counterbalance gravity, and the system collapses. This is 
the reason for the occurrence of density cusps in the CDM 
model for which P = 0. By contrast, a logotropic equa¬ 
tion of state can sustain gravity and develops flat density 
cores instead of cusps. The fact that the logotropic equa¬ 
tion of state ( |49| yields a Lane-Emden equation of index 
n = — 1 clearly demonstrates that it can be interpreted 
as the limit of a polytropic equation of state with 7 —> 0 
(i.e. n —> — 1) and K —> 00 such that K 7 = A is fi¬ 
nite |53j. In a sense, the logotropic equation of state is a 
“regularization” of the polytropic equation of state with 
index n = — 1. 


20 A self-gravitating system described by a distribution function of 
the form /(r, v) = /(e), where e = v 2 /2 + <3? is the individual en¬ 
ergy of the particles, has a barotropic equation of state P(p) ob¬ 
tained by eliminating 4>(r) from the relations p = f f dv = p(4>) 
and P = (1/3) f fv 2 dv = P(4>) j60]. For example, a system 
described by a polytropic distribution function has a polytropic 
equ ation of state m- In a sense, the logotropic equation of state 
( |49| ) is associated with a Lorentzian distribution function E3I- 
However, the connection is not direct because the Lorentzian is 
not normalizable in d = 3, except if we introduce a maximum 
bound on the velocity. Therefore, rigorously speaking, the lo¬ 
gotropic equation of state cannot be derived from a distribution 
function and, consequently, P/p does not represent a velocity dis¬ 
persion cr 2 . It is more likely that the logotropic equation of state 
comes from a field theory at T = 0 as discussed in Appendix [Cl 
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FIG. 17: Normalized density profile of a logotropic sphere. It 
decays as when f —> +oo. 


Since the logotropic spheres are homologous (53i, they 
generate a univeral density profile. Indeed, if we rescale 
the density by the central density po and the radial dis¬ 
tance by the core radius ?'o, we get an invariant profile 
l/0(£). The normalized density profile of a logotropic 
sphere is represented in Fig. [lT] The density of a lo¬ 
gotropic sphere decreases as at large distances oTSl . 
This can be compared to the r ~ 2 behavior of the isother¬ 
mal sphere d]. 21 Because of this slow decay, the total 
mass of a logotropic sphere is infinite. This implies that 
the logotropic equation of state cannot hold at large dis¬ 
tances. We note that the empirical Burkert density pro¬ 
file, that provides a good fit of many dark matter halos, 
decreases at large distances as r -3 .'T9 . This is sub¬ 
stantially steeper than the r _1 behavior of a logotropic 
halo. In practice, the finite mass and finite radius of dark 
matter halos result from tidal effects or from incomplete 
violent relaxation. Tidal effects and incomplete violent 
relaxation alter the density profile at large distances and 
steepen it. A famous example of distribution functions 
taking tidal effects into account is the King model jf78] . 
The King model at the limit of microcanonical stability 
can be approximated by the modified Hubble profile that 
has a core and that decreases at large distances as r -3 , 
similarly to the Burkert profile m- Since tidal effects, or 
other complicated effects such as incomplete relaxation, 
are not taken into account in the logotropic model, we 
should not give too much credit on the asymptotic behav¬ 
ior of its density profile at very large distances. 22 How¬ 
ever, it is relevant to mention that, in a recent paper, 


21 The singular logotropic sphere is p s = (A /8nG) 1 / 2 r 1 l53l and 
the singular isothermal sphere is p s = (l/2irG/3m)r - 2 [zzl. 

22 In the context of hierarchical clustering (small gravitationally 
bound clumps of dark matter halos merge to form progressively 
larger objects), we suggest that the core of the objects after merg¬ 
ing remains logotropic while their halo develops a density profile 
decreasing as r~ 3 due to tidal effects and incomplete relaxation. 


Burkert [53] observes that the slope of the density profile 
of dark matter halos close to the core radius is approxi¬ 
mately equal to —1 (see the upper right panel of his Fig. 
1) which is precisely the density exponent of a logotrope. 
This is a first hint that a logotropic equation of state may 
be relevant in the case of dark matter halos. 


Using the Lane-Emden equation (1281, the mass profile 


M(r) = f 0 p{r')4:TTr' dr 1 is given by 


M(r) = 47rp 0 r 3 £ 2 6l'(£). 


(130) 

The circular velocity defined by v 2 (r) = GM (r )/r can be 
expressed as 


Vc(r ) = ^TrGporfed'tt). 


(131) 


We define the halo radius ru as the radius at which 
p/po = 1/4- The dimensionless halo radius is the so¬ 
lution of the equation 9(£h) = 4. We numerically find 
£ h = 5.8458 and 0'(£ h ) = 0.69343. Then, r h = £ h r 0 . The 
normalized halo mass at the halo radius is given by 


M h 

Pod 


v 2 c (r h ) _ fl'(fr) 

Gpod a 


= 1.49. 


(132) 


This value is relatively close to the value Mh/pod = 1-60 
USED] obtained with the Burkert profile defined by 


P(r) 

Po 


(1 + a-)(l + x 2 ) ’ 


r 

x = —, 
Th 


(133) 


1.98 


ln(l + x) — arctana; + - ln(l + x 2 ) 


-I 1/2 


v c (r) 
v c (rh) 

(134) 

The density and circular velocity profiles of a logotrope 
are plotted as a function of the radial distance normal¬ 
ized by the halo radius in Figs. [18] and [19] They are 
compared to the observational Burkert profiles of Eqs. 


(1331 and (134). We see that the logotropic model gives 


a good agreement with the Burkert profile up to the halo 
radius, i.e. for r < rh- At larger distances, the agreement 
is less good because, in the logotropic model, the density 
decreases too slowly. This implies that the circular ve¬ 
locity keeps increasing after rh while the circular velocity 
obtained from the Burkert density profile reaches a max¬ 
imum and decreases at large distances. However, there 
exists galaxies such as LSB galaxies presenting rotation 
curves in which the circular velocity increases monoton- 
ically. On the other hand, as already indicated, the lo¬ 
gotropic model is not expected to be valid at large dis¬ 
tances since it does not take tidal effects or incomplete 
violent relaxation into account. Therefore, the fact that 
it can account for the observations up to the halo radius 
is already satisfying. This is sufficient to determine the 
physical characteristics of dark matter halos. 

The logotropic equation o f stat e has a very interesting 
property. According to Eq. (127), we have 


Porh 


= ( — ) 
\AttGJ 


1/2 


€h- 


(135) 
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r/r, 

n 


FIG. 18: Normalized density profile of a logotropic sphere 
(solid line) compared with the Burkert profile (dashed line). 



FIG. 19: Circular velocity profile of a logotropic sphere (solid 
line) compared with the Burkert profile (dashed line). 


Therefore, if the logotropic temperature A is the same 
for all the halos (we claim in Sec. VI that A is a univer¬ 
sal constant), the surface density So = po^h is also the 
same. This is precisely what is observed. Indeed, it is 
an empirical fact that the surface density So = porh is 
approximately the same for all the galaxies even if their 
sizes and masses vary by several orders of magnitude (up 
to 14 orders of magnitude in luminosity). Its best-fit 
value is E 0 = Af 0 /pc 2 [S&HSE]- To our knowledge, 

this observation has not been explained theoretically. It 
turns out that this result is a direct consequence of the 
logotropic equation of state if we assume that the dark 
matter halos have the same logotropic constant A. This 
is particularly appealing if we interpret this constant as 
a temperature in a generalized thermodynamical frame¬ 
work (see Appendix C5c). In this interpretation, the 
universe is “isothermal” except that isothermality does 
not refer to a linear equation of state but to a logotropic 
equation of state. 


The relation (135) allows us to determine the lo¬ 


gotropic temperature A from the measurement of the 
surface density E 0 . We find 


A = 4t tG ( = 2.13 x 1CT 9 


gm 1 s 2 . 


(136) 


The dimensionless logotropic temperature is then given 
by 


B = 


4t tG 


£q 1 — 12 


m, 0 


) = 3.53 x 10" 3 . (137) 


We note that this value satisfies the cosmological bound 
0 < B < 0.09425 obtained in Sec. |IV F[ as well as the 
bounds 0 < B < 0.0262 and 0 < B < 0.0379 obtained 
from the measurement of the shift parameter from CMB 
in Sec. |IVM| This agreement is particularly rewarding 


because the value of B obtained in Eq. (|137|) is based on 
galactic observations that are complete 
the cosmological observations of Sec. 


y different from 


IV F and IV M 


Therefore, a logotropic equation of state can simultane¬ 
ously describe dark matter halos and account for cosmo¬ 
logical constraints. This may be a hint that dark matter 
and dark energy are the manifestations of a single dark 
fluid. 

There are interesting consequences of this result. Ac¬ 
cording to Eq. (132), the mass of the halos calculated at 
the halo radius is given by M h = 1.49E 0 r 2 . Since the 
surfa ce density of the dark matter halos is constant, Eq. 


(132) implies that M/,/Af 0 = 210(r/ l /pc) 2 oc r 2 . This 
scaling is consistent with the observations [80]. On the 
other hand, Strigari et al. [55 have proposed that all 
dwarf spheroidal galaxies (dSphs) of the Milky Way have 
the same total dark matter mass contained within a ra¬ 
dius of r u = 300 pc. From the observations, they obtain 
log(M3oo/Af 0 ) = 7.0 ^q |. This result is still unexplained. 
We now show how the logotropic equation of state im¬ 
mediately leads to this result. According to Eq. (130), 
using £ u = r u /r 0 , we have 


(138) 


M 300 = 4irpor 0 rl6' ( ^ 


Introducing the halo radius Vh = o, we obtain 


M 3 oo = ±Kpor h -±# ( th- 
\ r h 


(139) 


As we have already indicated, it is an observational ev¬ 
idence that the surface density pcT/i of the ha los is a 
constant. 23 In principle, the last term in Eq. (139) is 


not a constant since it depends on rh which substantially 
changes from halo to halo. However, for the logotropic 


23 In our model, this c an b e explained theoretically by the univer¬ 
sality of A (see Sec. VI i. 





























22 


distribution, we have the asymptotic result 0(£) ~ (,/\/2 
for £ —► +oo [53]. More precisely, for £ > 8 (typically), 
0'(£) undergoes damped oscillations about l/-\/2. There¬ 
fore, 6'(£) quickly reaches a constant value 1 /v/2- For the 
dSphs considere d in .ESI , £,hr u /r h > 1 (see, e.g., Table 2 
of [81]), so Eq. (139) may be replaced by 


M 30 o = 


471-Eq rg 
£/»\/2 


(140) 


which is a constant in agreement with the claim of Stri- 
gari et al. [52!. Furthermore, using E 0 = 141M 0 /pc 2 , 
the numerical application gives 


M 3nn = 1.93 x 10 7 M 


0i 


(141) 


leading to log(M 3O o/M 0 ) = 7.28 in very good agreement 
with the observational value. It is also relevant to express 
this result directly in terms of the logotropic temperature 
A under the form 


M 300 


1/2 


(142) 


which shows that the constancy of M 30 o is due to the 
universality of A. 


VI. A PREDICTION OF THE LOGOTROPIC 
TEMPERATURE 

In the previous section, the values of A and B have 
been deduced from the observations. In this section, we 
show that they can actually be predicted from general 
considerations. 

Introducing the cosmological density e A = p A c 2 = (1 — 
fI m> o)eoj which corresponds to the present value of the 
dark energy term (internal energy) in Eq. (571, we can 
rewrite Eq. (601 as 


P*C _ f^m.O jB 

CA 1 0 


(143) 


Since B is expected to be very small, this ratio is huge. 
This makes us think of the huge ratio pp/p A ~ 10 123 be¬ 
tween the vacuum energy, associated to the Planck den¬ 
sity pp = 5.16 x 10" gm -3 , and the dark energy (or 
cosmological constant A), associated to the cosmological 
density p A = 6.72 x 10 -24 gm -3 . In the framework of 
the ACDM model, this huge ratio is at the origin of the 
cosmological constant problem [241 25]. In the present 
approach, in which there is no cosmological constant nor 
dark energy, it has a completely different significance. We 
propose that it determines the dimensionless logotropic 
temperature B. Therefore, we propose to identify the 
density p* in the logotropic equation of s tate ( |49| ) with 
the Planck density pp. In that case, Eq. (143) becomes 


Pp 

PA 


f^m.O 1+1/B 


1 - a 


m, 0 


(144) 


This yields approximately B ~ 1/ \n(pp/p A ) ~ 

1/[123ln(10)], where 123 is the famous number occur¬ 
ring in the ratio Pp/pa ~ 10 123 . More precisely, 


B = 


1 


In (- 1 

\ “m,o Pa ) 


= 3.53 x 10 


-3 


(145) 


This value turns out to be in perfect agreement with 
the value obtained from the measurement of the surface 
density of dark matter halos in Sec. [Vj If we accept 
that this agreement is not fortuitous, or coincidental, this 
is a strong argument in favor of the logotropic model. 
Therefore, by taking p* = pp in the logotropic equation 
of state (49), we obtain a prediction of the parameter 


B which is in perfect agreement with the observations. 
With this value, there is no free parameter in our model, 
so we can predict all the measurable quantities that occur 
in cosmology a nd in the study of dark matter halos (see 
Table I in Sec. VII). We emphasize that our estimate of 


B only depends on the values of pp , p\ and that 

are known accurately from the observations. 


In conclusion, the logotropic equation of state (49) can 
be written as 


P = Bca In ( — 
, Pp 


(146) 


with pp = 5.16 x 10" g m -3 
gm -1 s" 
temperature is 


10 7 gm 1 s 2 and B = 3.53 x 10 3 


6a — (1 — fi m ,o) e o — 6.04 x 
The logotropic 


A = Be A = 2.13 x 10 -9 gm -1 s -2 . 


(147) 


It is of the order of the cosmological density e A divided 
by ln(p P /p A ) ^ [123ln(10)]. 

The ACDM model is recovered for B = 0. This cor¬ 
responds to pp -a + 00 , hence h —► 0. Therefore, the 
fact that B is small but nonzero shows that quantum 
mechanics plays a role in the late universe in relation to 
dark energy. 


The logotropic equation of state (1461 has several nice 
properties. 

(i) If our approach is correct, there is no dark mat¬ 
ter nor dark energy. There is just one dark fluid de¬ 
scribed by the equation of state (|146|. In that case, 


it may be relevant to interpret the logotropic temper¬ 
ature A = 2.13 x 10~ 9 gm _1 s^ 2 as a fundamental con¬ 
stant which supersedes the cosmological constant. We 
note that it depends on all the fundamental constants of 


physics h, G, c, and A [see Eqs. (145) and (147)]. Then, 
Eq. (56) with p* = pp may be used to predict fI mi o for a 
given value of the present energy density eo determined 
from the observations by the Hubble constant. In the 
present interpretation fl m ,o does not represent the pro¬ 
portion of dark matter. It simply gives the proportion 
of the rest-mass energy of the dark fluid as compared 
to its total energy (the remaining energy being internal 
energy). 
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(ii) Since dark matter and dark energy are the mani¬ 
festation of a single dark fluid, there is no cosmic coin¬ 
cidence problem. On the other hand, the cosmological 
constant probl em p p/pa ~ 10 123 is translated into an 
equation [Eq. (144)] that determines the logotropic tem¬ 
perature B. 


(iii) In the polytropic model (see Sec. IIHi, the pres¬ 


sure is negative because K < 0, i.e. because the poly¬ 
tropic temperature is negative. In the logotropic model, 
the logotropic temperature A is positive, and the pres¬ 
sure is negative because p < pp. Positive temperatures 
are more satisfying on a physical point of view than neg¬ 
ative temperatures. On the other hand, while the pres¬ 
sure P(p) is negative, which is required at the cosmo¬ 
logical scale to produce an acceleration of the expansion 
of the universe, its derivative P'(p) is positive, which is 
required at the scale of dark matter halos to counteract 
self-gravity and avoid density cusps. 

(iv) It is interesting to see that both the cosmologi¬ 
cal density p\ and the Planc k de nsity pp appear in the 
logotropic equation of state (146). The presence of the 


Planck density is rather unexpected, and intriguing, be¬ 
cause the logotropic equation of state only describes a 
relatively old universe, dominated by dark matter and 
dark energy (pa) well after the phase of early inflation 
dominated by quantum mechanics (pp). The fact that 
the Planck den sity appears as a density scale in the log¬ 
arithm of Eq. (1461 is very interesting because only the 


logarithm of the ratio of the cosmological density over 
the Planck density has a physical meaning since these 
quantities differ by 123 orders of magnitude. The occur¬ 
rence of the Plan ck d ensity suggests that the logotropic 
equation of state (1461 may be the limit of a more general 


equation of state connecting the very early universe (pp) 
to the very late universe (pa)- In relation to this obser¬ 
vation, we note that the pressure in the logotropic model 
vanishes precisely at the Planck scale, so our treatment 
breaks down at that scale. 

Now that B and A have been derived from theoret¬ 


ical considerations [see Eqs. (145) and (147)], we can 
reverse the arguments of S ec. VI and predict the values 


of E 0 and M 300 [see Eqs. (1351 and (142)]. We obtain 
E 0 = 141M 0 /pc 2 and M 300 = 1.93 x 10 7 M 0 in per¬ 
fect agreement with the observations, and without any 
free parameter. We can also obtain an estimate of the 
Jeans length A j at the beginning of the matter era where 
perturbations start to grow. We assume that the matter 
era starts at Ri = 10 -4 , corresponding to the epoch of 
matter-radiation equality. In this era, we can make the 
approximation e = pc 2 , so the Jeans wavenumber is given 
by m- 


k 2 j = 


4ttG P R 2 


(148) 


From Eqs. (66) and (97), we find pi = 2.54 x 10 _12 g/m 3 
and (c 2 /c 2 )i = 9.33 x 10~ 15 . This leads to a Jeans length 
A j = 2Tr/kj = 1.25 x 10 18 m = 40.4pc which is of the 
order of magnitude of the smallest known dark matter 


halos such as Willman I (th = 33pc), see, e.g., Table 2 
of |8Tj. We predict that there should not exist halos of 
smaller size since the perturbations are stable for A < A j. 
This is in agreement with the observations. By contrast, 
in the CDM model, since P = 0, the velocity of sound 
c s =0. As a result, the Jeans length is zero (Aj = 0), 
implying that the homogeneous background is unstable 
at all scales so that halos of any size should be observed 
in principle, which is not the case. Therefore, a small but 
nonzero value of B , yielding a nonzero velocity of sound 
and a nonzero Jeans length, is able to account for the 
minimum observed size of dark matter halos. It also puts 
a cut-off in the density power spectrum of the perturba¬ 
tions and sharply suppresses small-scale linear power. Fi¬ 
nally, it solves the cusp problem and the missing satellite 
problem. On the other hand, if we make the numerical 
application for the present-day universe (Rq = 1), we find 
po = 2.54 x 10~ 24 g/m 3 , (c 2 /c 2 )o = 9.42 x 10~ 3 and A j = 
1.25 x 10 26 m = 4.06 10 9 pc. The Jeans length is so large 
(of the order of the horizon Xp = c/Hq = 1.32 x 10 26 m 
since A j ~ c/y/Gpo ~ c/Hq ) that there is no Jeans in¬ 
stability in the present universe. A detailed study of the 
perturbations in the logotropic model is beyond the scope 
of this paper and will be considered in a future work. 


VII. NUMERICAL APPLICATIONS 


In this section, we provide the values of the princi¬ 
pal quantities that occur in cosmology and in the study 
of DM halos (see Table I). We make the numerical ap¬ 
plication for: (i) B = 0, corresponding to the ACDM 
model, (ii) B = 0.0942 5, cor responding to the cosmolog¬ 
ical constraint of Sec. IVF and (iii) B = 3.53 x 10 -3 , 
corresponding to the prediction of Sec. [Vi] 

We recall that (R w , e w ,t w ) refer to the values at which 
the pressure of the universe becomes negative, (i? c , e c , t c ) 
refer to the values at which the universe accelerates, (I? 2 , 
e 2 , f 2 ) refer to the transition between dark matter and 
dark energy (i.e. when the internal energy of the dark 
fluid dominates its rest-mass energy), to is the age of 
the universe, q 3 is the present value of the deceleration 
parameter, w 3 is the present value of the equation of 
state parameter, c Sy o is the present velocity of sound, 
(R a , e. s , t s ) refer to the values at which the velocity of 
sound exceeds the speed of light, ( Rm , cm, <m) refer 
to the values at which the universe becomes phantom, 
R is the shift parameter from CMB, A is the logotropic 
temperature, £o is the surface density of DM halos, and 
M 3 qo is the mass of dwarf halos enclosed within a sphere 
of radius r u = 300 pc. 

We note that for the predicted value of B = 3.53 x 
10~ 3 , the cosmological parameters are almost the same 
as in the ACDM model. This is satisfactory since the 
ACDM model works well at the cosmological scale. In 
this sense, there is no important difference, from the 
observational viewpoint, between the logotropic model 
(B = 3.53 x 10~ 3 ) and the ACDM model (B = 0) at 
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B = 0 (ACDM) 

B = 0.09425 

B = 3.53 x 10~ 3 

Rw 


2.09 x 10~ 2 

7.00 x 10" 42 

(e/eo)™ 


3.02 x 10 4 

7.97 x 10 122 

t w (Gyrs) 


0.0535 

3.29 x 10" 61 

Rc 

0.574 

0.576 

0.574 

(e/eo)c 

2.17 

2.05 

2.17 

t c (Gyrs) 

7.18 

7.33 

7.19 

R 2 

0.723 

0.744 

0.723 

( £ / e 0)2 

1.45 

1.33 

1.45 

t 2 (Gyrs) 

9.61 

10.1 

9.63 

to (Gyrs) 

13.8 

14.0 

13.8 

4o 

-0.589 

-0.692 

-0.593 

w 0 

-0.726 

-0.794 

-0.729 

(c 2 s/c 2 ) 0 

0 

0.331 

9.42 10 -3 

Rs 


1.26 

3.77 

(e/eo )s 


0.910 

0.741 

ts (Gyrs) 


17.3 

34.5 

Rm 


1.59 

4.75 

(e/eo )m 


0.889 

0.7405 

t M (Gyrs) 


20.7 

38.3 

R 

1.7342 

1.7505 

1.7348 

A (g m -1 s -2 ) 


5.69 x 10~ 8 

2.13 x 10 -9 

E 0 (M 0 /pc 2 ) 


1523 

141 

M 300 (Mq) 


2.08 x 10 8 

1.93 x 10 7 


TABLE I: Quantities that occur in cosmology and in the study of dark matter halos for different values of B (see the text for 
details). Some observational values are R c = 0.571 ± 0.013 [73], qo = — 0.53lg'jg jS2], and R 2 = 0.719 ± 0.017 [73). 


the cosmological scale. Differences will appear in the re¬ 
mote future, in about 20 Gyrs, when the velocity of sound 
approaches the speed of light and the universe tends to 
become phantom (the logotropic model may break down 
before). However, on a theoretical point of view, the 
logotropic model solves the cosmic coincidence problem 
(since there is just one fluid) and avoids the cosmologi¬ 
cal constant problem (it gives an interpretation to the 
number 123 in terms of a logotropic temperature B ). 
On the other hand, a small but nonzero value of the 
logotropic temperature B = 3.53 x 10~ 3 is very impor¬ 
tant at the scale of dark matter halos. First, it yields a 
nonzero velocity of sound and a nonzero Jeans length that 
solve the missing satellite problem and the cusp prob¬ 
lem. Furthermore, the value of the Jeans length at the 
beginning of the matter-dominated era is consistent with 
the observed size of the smallest known halos, and ex¬ 
plains therefore the existence of a minimum halo radius 
Rmin ~ 10 pc in the universe. On the other hand, the 
value B = 3.53 x 10~ 3 of the logotropic temperature can 
account remarkably well for the structural properties of 
dark matter halos. Not only our model explains why the 
surface density E 0 and the mass M 300 of DM halos are 
the same for all the halos but it also predicts their values 
with high accuracy (see Table I). Furthermore, it shows 
that their values are related, through A and B , to the 


Planck de nsity pp an d to t he cos molo gical density p\ 
[see Eqs. (1351, (142), (145), and (147)]. Such a depen¬ 
dence was unexpected and remains relatively mysterious. 
Our model first suggests that dark matter and dark en¬ 
ergy are the manifestation of a unique entity (dark fluid) 
since the properties of dark matter halos such as £o and 
M 30 o are related, through the logotropic temperature A, 
to cosmological observables such as p\ (dark energy). On 
the other hand, the fact that the Planck density enters 
in the logotropic equation of state initially designed to 
model dark matter and dark energy is intriguing. It sug¬ 
gests that quantum mechanics manifests itself at the cos¬ 
mological scale in relation to dark energy. This may be a 
hint for a fundamental theory of quantum gravity. This 
also suggests that the logotropic equation of state may 
be the limit of a more general equation of state providing 
a possible unification of dark energy (p\) in the late uni¬ 
verse and inflation (vacuum energy pp) in the primordial 
universe. This may be a clue for the elaboration of a 
more general theory incorporating inflation. 


VIII. CONCLUSION 

We have proposed that the universe is made of a sin¬ 
gle dark fluid described by a logotropic equation of state. 
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In our model, there is no cosmological constant, no dark 
matter, and no dark energy. As a result, the cosmological 
constant problem and the cosmic coincidence problem do 
not arise. What we traditionally call dark matter, dark 
energy, and cosmological constant can be related to the 
intrinsic properties of the dark fluid. Dark matter corre¬ 
sponds to the rest-mass energy of the dark fluid and dark 
energy corresponds to its internal energy. The logotropic 
temperature A of the dark fluid can be interpreted as a 
fundamental constant of physics superseding the cosmo¬ 
logical constant. 

We have first determined bounds on the dimension¬ 
less logotropic temperature B by using cosmological con¬ 
straints: 

(i) The condition that the velocity of sound is not al¬ 
ways imaginary implies B > 0. 

(ii) The condition that the present velocity of sound 
is real implies B < B\ = 0.377. This condition coin¬ 
cides with the condition that the present universe is non¬ 
phantom. 

(iii) The condition that the present velocity of sound is 
less than the speed of light implies B < B±/2 = 0.1885. 

(iv) The condition that the present universe is non- 
relativistic implies B < i?i/4 = 0.09425. 

(v) The measurement of the shift parameter from CMB 
implies B < 0.0379 or B < 0.0262 depending on the 
estimates. 

We have then proposed that dark matter halos are de¬ 
scribed by the same logotropic equation of state. For 
B > 0, this equation of state leads to flat cores instead 
of density cusps, so it solves the cusp-core problem. Fur¬ 
thermore, it generates a universal density profile and a 
universal rotation curve that agree with the observational 
Burkert profiles (32 U P to the halo radius 77. 24 More 
specifically, and more strikingly, a logotropic equation 
of state explains naturally (i) the recent observation of 
Burkert [55] that the density of dark matter halos de¬ 
creases as r -1 close to the core radius, (ii) the observation 
that the surface density of dark matter halos is the same 
for all the halos [5HH55] , and (iii) the observation that the 
mass of dark matter halos contained in a fixed radius (e.g. 
300 pc) is the same for all the dwarf halos [59] . Using the 
measured value of the surface density E 0 = 141M 0 /pc 2 
I55H55] . we obtain a value B = 3.53 x 10 -3 for the lo¬ 
gotropic temperature. From this value, we find that the 
mass of dark matter halos contained in a radius of 300 
pc is M 30 o = 1.93 x 10 7 M 0 , in agreement with the value 
obtained by Strigari et al. j5Uj - 

Finally, we have argued that the reference density p* 


24 The rotation curve of logotropic dark matter halos continues to 
increase (instead of decreasing) at larger distances. This is con¬ 
sistent with the rotation curves of certain galaxies, such as low 
surface brightness (LSB) galaxies, that are particularly well iso¬ 
lated. For other galaxies, tidal effects or complex physical pro¬ 
cesses (e.g. incomplete relaxation) have to be taken into account 
at large distances to confine the system. 


appearing in the logotropic equation of state should be 
identified with the Planck density pp. This predicts the 
value of the logotropic temperature to be B = 3.53 x 10 ~ 3 
in perfect agreement with the value deduced from the 
observations. The corresponding value of the dimen¬ 
sional logotropic temperature is A = Be a = 2.13 x 
10~ 9 gm - 1 s -2 . Such a small value of B yields sensi¬ 
bly the same results as the ACDM model (B = 0) at the 
cosmological scale. However, having B = 3.53 x 10 ~ 3 
instead of B = 0 is crucial at the scale of dark matter 
halos to explain their properties, as shown above. Fur¬ 
thermore, it implies a nonzero velocity of sound and a 
nonzero Jeans length that is, at the beginning of the mat¬ 
ter era, of the order of the minimum size (~ 10 pc) of the 
observed dark matter halos. This also solves the missing 
satellite problem. 

The next step is to justify the logotropic equation of 
state from first principles. We note that our approach is 
already very economical. Instead of having dark matter 
and dark energy with two different equations of state, 
we just have one dark fluid with one equation of state. 
The logotropic equation of state involves two parameters, 
the reference density p* and the logotropic temperature 
A. They are related to each other by Eq. (56) which 
depends on the known (observed) values of eo and ft m 0 . 
Assuming that p* = pp, this equation determines A, so 
there is no free parameter in our model. The polytropic 
equation of state (Sec. |IIH ) also involves two parameters, 
the polytropic index 7 and the polytropic temperature K. 
They are related to each other by an equation similar to 
Eq. (56). However, it seems difficult to justify why the 
polytropic index should have a particular value such as 
7 = —0.089 [35]. On the other hand, in the logotropic 
model, the logotropic temperature A is positive so it can 
be related to an energy scale A = 3.53 x 10 _ 3 eA- In 
the polytropic model, the polytropic temperature K is 
negative so its physical interpretation is not direct. 

We have proposed different possibilities to justify the 
logotropic equation of state from first principles: 

(i) The logotropic equation of state is related to Tsal- 
lis generalized thermodynamics ]54| . It is associated with 
the Log-entropy [53] (see Appendix C 5 c) which is a sort 
of regularized Tsallis entropy with index q = 0 (i.e. 7 = 0 
in the polytropic terminology). This marginal index, 
where power laws degenerate into a logarithm, may be 
viewed as a fixed point, or as a critical point, in the frame¬ 
work of Tsallis generalized thermodynamics. Therefore, 
if our cosmological model is correct, it would be a nice 
confirmation of the interest of generalized thermodynam¬ 
ics in physics and astrophysics. 

(ii) We have proposed that the dark fluid could be a rel¬ 
ativistic self-interacting scalar field representing, for ex¬ 
ample, a self-interacting BEC at T = 0. In that context, 
the logotropic equation of state arises from the GPequa- 
tion with an inverted quadratic potential [Eq. (C56)], or 
from the KG equation with a logarithmic potential [Eq. 

(C58)]. These potentials may be simpler to justify from 
fundamental physics than power-law potentials with ex- 
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ponents 2(7 - 1) = -2.18 [Eq. fC42| ] and 2 7 = -0.178 
[Eq. (C44 )1 associated with a polytropic equation of state 
with-T^- 0.089 BS|. 


If the equation of state is prescribed under the form 


(iii) In Sec. Ill we have mentioned that our ap¬ 


proach provides a new justification of the Cardassian 
model [68j. Inversely, the original justification of the 
Cardassian model, namely that the term zz(p) arising 
in the modified Friedmann equation (44) may result 
from the existence of extra-dimensions, could also be 
a way to justify the logotropic model corresponding to 
v(p) = ~(8 ttGA/3c 2 ) [In (p/pp) + 1]. In Appendix | b| we 
present a simple argument to justify the logotropic equa¬ 
tion of state and we relate this equation of state to a 
long-range confining force that could be a fifth force or 
an effective description of higher dimensional physics. 

In conclusion, the logotropic equation of state may be 
a good candidate for the unification of dark matter and 
dark energy. Indeed, it provides a good description of the 
cosmological evolution of the universe as a whole and, at 
the same time, accounts for many properties of dark mat¬ 
ter halos, some of them being until now unexplained. The 
logotropic equation of state is rather unique in account¬ 
ing for all these observational results in a unified manner. 
It is related to generalized thermodynamics and corre¬ 
sponds to the Log-entropy which is a sort of degenerate 
Tsallis entropy with index q = 0. It is also related to the 
Cardassian model. Finally, it can be given a field the¬ 
ory interpretation as it arises from a nice simplified form 
of GP and KG equations corresponding to an inverted 
quadratic potential or to a logarithmic potential, respec¬ 
tively. All these elements suggest that dark matter and 
dark energy may be the manifestation of a unique dark 
fluid with a logotropic equation of state. Some develop¬ 
ments of this model will be given in future works. 


Appendix A: Relativistic thermodynamics 

The local form of the first law of thermodynamics can 
be expressed as 


d I - ) = -Pd ( - ) + Td [ — 


(3 


(Al) 


where p = nm is the mass density, n is the number den¬ 
sity, and s is the entropy density in the rest frame. For a 
system at T = 0, or for an adiabatic evolution such that 
d(s/p) = 0, 25 the first law of thermodynamics reduces to 

P 


de = 


-dp. 


(A2) 


For a given equation of state, Eq. (A2) can be integrated 


to obtain the relation between the energy density e and 
the rest-mass density. 


25 The evolution of a perfect fluid is adiabatic. In particular, it can 
be explicitly checked that the Friedmann equations conserve the 
entropy m- 


P = .P(e), Eq. (A2) can be immediately integrated into 

de 


In p = 


P(e) 


(A3) 


If, as an example, we consider the “gamma law” equation 

of state [531 Ell : 


P = ( 7 -l)e 


we get 


P = Kp\ 


€ = 


K 7 

- TP ’ 

7 -l 


(A4) 


(A5) 


where K is a constant of integration. 

We now assume that the equation of state is pr escribed 
under the form P = P(p). In that case, Eq. (A2) reduces 
to the first order linear differential equation 


de 1 P{p) 

dp p p 


(A6) 


Using the method of the variation of the constant, we 
obtain 


e = Ape 2 + p 




(A7) 


where A is a constant of integration. 

As an example, we consider the polytropic equation of 
state 


7=1 + —. 
n 


P = Kp\ 

For 7 = 1, we get 

e = Ape 2 + Kp In p. 

For 7 7^ 1, we obtain 

K 

e = Ape 2 H- p 1 = Ape 2 + nP. 

7-1 


(AS) 

(A9) 

(A10) 


Taking A = 0, we recover Eqs. (A4) and (A5). Taking 


A = 1, a choice that we shall make in the following, we 
obtain Eqs. (38) and (39). 


We first assume n > 0 (i.e. 7 > 1). For p —>■ 0, 

e - pc\ P ~ K{e/e 2 )\ (All) 

and for p —» +00, 

e ~ nKp 1 , P ~ e/n ~ (7 — l)e. (A12) 

We now assume n < 0 (i.e. 7 < 1). For p —> 0, 

e ~ nKp 1 , P ~ e/n ~ (7 — l)e, (A13) 

and for p —» +00, 

e ~ pc 2 , P ~ K{e/c 2 y. (A14) 
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For a general equation of state P(p), we obtain 


e = pc + p 


P(p') 


dp' = pc 2 + u(p), 


(A15) 


where the primitive is determined such that u(p) does 
not contain terms in pc 2 . We note that u(p) may be 
interpreted as an internal energy density (see [851 and 
Appendix |C 4h . Therefore, the energy density e is the 


sum of the rest-mass energy pc 2 and the internal energy 
u(p). The rest-mass energy is positive while the internal 
energy can be positive or negative. Of course, the total 
energy e = pc 2 + u(p) is alwa ys po sitive. 

Remark: according to Eq. (A2), the velocity of sound 


defined by c 2 /c 2 = P'(e) satisfies the identity 


5 = ^ 

r 2 de 
dp 


(A16) 


Appendix B: A simple argument to justify the 
logotropic equation of state and its relation to a fifth 

force 


We have seen in Sec. |II H| that the energy density e of 
a relativistic dark fluid at T = 0 (or a perfect fluid) is 
the sum of two terms: a rest-mass term pc 2 mimicking 
dark matter and an internal energy term 


u(p) = P J 


P(p') , , 

—yrdp 


(Bl) 


mimicking a “new fluid”. If we want the new fluid to 
mimick dark energy, we should have P = —u. Substitut¬ 
ing the relation u(p) = —P{p) in Eq. (Bl) and solving 
the resulting differential equation, we find P = cst. which 


returns the model of Sec. IIF equivalent to the ACDM 


model. However, we have explained that a constant pres¬ 
sure does not account for the structure of dark matter ha¬ 
los. Therefore, a next guess is P = —u — A where A is a 
const ant . Substituting the relation u(p) = —P(p) — A in 
Eq. (Bl) and solving the resulting differential equation, 


we find the logotropic equation of state P = Aln(p/p*). 
We note that P ~ — u as p —> 0, similarly to the equa¬ 
tion of state of dark energy. Another possible guess 
is P = u/n with n —> —1. Substituting the relation 
u(p) = nP(p) in Eq. (Bl) and solving the resulting 


differential equation, we find the polytropic equation of 
state P = Kp 1+1 t n . For n = —1, we obtain P = K. 
The polytropic equation of state with an arbitrary index 
n is interesting at a general level (beyond dark energy). 
For example, the index n = 1, leading to a quadratic 
equation of state P = Kp 2 , corresponds to a stiff fluid 
(P = it) possibly existing in the very early universe [44]. 

The pressure P can be associated with a self¬ 
interaction potential between particles. If we consider 
a power-law interaction of the form 


U(nj) = U 0 r° 


(B2) 


corresponding to a force F{r l: j) = —aUor°/ 1 , the virial 
theorem writes (see, e.g., Appendix I of [55]'): 

2 K-aW = 3 PV, (B3) 

where K is the kinetic energy, W is the potential energy, 
P is the pressure, and V is the volume of the system. 
Taking K = 0 for a dark fluid at T = 0 and identifying 
W/V as the internal energy u, we obtain (see an alterna¬ 
tive derivation of this result in the Appendix of [65]): 


P = —era. 
3 


(B4) 


The equation of state P = —u of dark energy corresponds 
to a = 3 yielding U = U 0 rfj and F = —3U 0 r 2 j. This is 
a long-range confining force that could be a fifth force 
or an effective description of higher dimensional physics 
[ 68 ] . The equation of state P = u/n associated with a 
polytrope of index n corresponds to a = —3/n yielding 
U = U 0 r~^ n and F = (3/n)U 0 r~j 3 ^ n ~ 1 . In particular, 
the equation of state of radiation P = u/3 associated 
with a polytrope of index n = 3 corresponds to a = 

— 1 yielding U = U 0 r/ '- 1 and F = U 0 r~- 2 (Coulombian 
force). The stiff equation of state P = u associated with a 
polytrope of index n = 1 corresponds to a = — 3 yielding 
U = U 0 r~ 3 and F = 3U 0 r~j 4 (attractive force ). 26 We 
note that the dark energy equation of state and the stiff 
equation of state correspond to a polytropic index n = 

— 1 and n = +1 respectively, and to a potential with an 
exponent +3 and —3. This symm etry may be related to 
the comment following Eq. (C56). 


For a logarithmic potential of interaction of the form 
U(rij) = U 0 \nrij, (B5) 

the virial theorem writes (see, e.g., Appendix I of [ 86 ]): 


2 K- 


UeN 2 


= 3 PV, 


(B 6 ) 


where N is the number of particles in the volume V. 
Therefore, it is tempting to associate the constant A 
(logotropic temperature) in the equation of state P = 
—u — A as arising from a logarithmic potential of inter¬ 
action of the form of Eq. (B5) such that 


A = 


U 0 N 2 

6V 


(B7) 


The total potential of interaction is therefore U = Vor/j + 
U 0 In nj and the total force is F = —3Uqt 2 a — Uq jr^. The 


meaning of Eq. (B7) is not clear, but it is interesting 


to note that A > 0 for an attractive potential Uq > 0 
(expected for dark energy), which is the correct sign of 
the logotropic temperature. 


26 Actually, in Zel’dovich model ED, the stiff equation of state 
arises from a screened Coulombian potential (see Sec. 4.2. of 

El). 
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Appendix C: Field theory 

1. The Klein-Gordon equation 

The Klein-Gordon (KG) equation for a complex scalar 
field <p writes 


1 d 2( t> _ 

c 2 dt 2 * 


h 2 


, 24- 

H- 2 

c z 


,+2 ^W (t> ~ 0 ' (C1) 


where V = V{\<p\ 2 ) is the self-interaction potential and <f> 
is an external potential that, in a simplified model, can be 
identified with the gravitational potential (see [581 for a 
fully general relativistic treatment). The energy density 
and the pressure of the scalar field are given by 


1 

2d 2 


P = 


2 c 2 


dcp 

~dt 


dcp 

at 


^ |v ^>| 2 + v(\4>\ 2 ) ■ 


2 h 2 


>r, (C 2 ) 


2 2 
m z c z 


Im'-vm*)- 2na 


|^| 2 . (C3) 


Mathematically, we can always make this change of vari¬ 
ables. However, we emphasize that it is only in the non- 
relativistic limit c —¥ +00 that ip has the interpretation of 
a wave function, and that \ip\ 2 = p has the interpretation 
of a rest-mass density. In the relativistic regime, ip and 
p = \tp\ 2 do not have a clear physical interpretation. We 
will call them “pseudo wave function” and “pseudo rest- 
mass density”. Nevertheless, it is perfectly legitimate to 
work with these variables [ 55 ;. 

Substituting Eq. (C8) in the KG equation (Cl), we 
obtain 


h 2 


dV 


d 2 ip dip h 2 . , , , 

wae“ ,ft a-ta A,A+m ® ,>+m 5pp 


ip = 0. (C9) 


On the other hand, the energy density and the pressure 
can be written in terms of ip as 


h 2 


2 m 2 c 2 


dip 


dt 


■ ,2 2 ft t ( 1 3^* 

+ Me 2 + —Im ip—- 
m \ ot 


h2 \Vip\ 2 + V(\ip\ 2 ), (CIO) 


2 m 2 


The KG equation without self-interaction can be viewed 
as the relativistic generalization of the Schrodinger equa¬ 
tion. Similarly, the KG equation with a self-interaction 
can be viewed as the relativistic generalization of the GP 
equation. In order to recover the Schrodinger and GP 
equations in the nonrelativistic limit c —> + 00 , we make 
the transformation 


<p{r, t) = Ae- imcH / h ip(Y,t), 


(C4) 


where A is a constant. This constant can be determined 


by the following argument. Substituting Eq. (C4) in the 
expression of the energy density (C2), we get 


A 2 ^ 

2d 2 


dip 

Ik 


m~c 

h 2 


A 2 m r ( , dip* 

-JT lm l^-dt 


-H 2 

+4 : IW’| 2 + k(|V'I 2 ). (C5) 


In the nonrelativistic limit c —> + 00 , we have e 
where p is the rest-mass density. On the other hand, 
according to Eq. (C5), 


e 


h 2 


■M 2 - 


(C6) 


If we interpret ip as the wavefunction normalized such 
that \ip\ 2 = p, we find by identification that 


A=l. 

m 


Therefore, Eq. (C4) can be rewritten as 

0(r,t) = -e~ lmc2t/h iP(r,t). 

m 


(C7) 


(C8) 


P = 


h 2 


2 m 2 c 2 


dip ft , / , d'ip* 

— + —Im ip—- 

ot nn \ ot 


h2 \S7ip\ 2 - V(\ip\ 2 ). 


2 m 2 


(Cll) 


2. The Gross-Pitaevskii equation 

Taking the nonrelativistic limit c —> +00 of the KG 
equation (C 91 , we obtain the nonlinear Schrodinger equa¬ 
tion 


dip h 2 dV 

in— = -—Aip + m9ip + m——^ip. 
dt 2 to dl^l 2 

It can be written as a GP equation of the form 


■ h dip h 2 


■ to 4-'0 4- mhd^l 2 )^ 


(C12) 


(C13) 


P c with a potential 


MIV’I 2 ) = 


dV 
d\ip\ 2 


i.e. h{p) = V\p). (C14) 


The GP equation describes a BEC at T = 0. The 
GP equation with a cubic nonlinearity (corresponding to 
h(p) oc p) can be derived from the mean held Schrodinger 
equation with a pair contact potential (see Sec. II.A. of 
[89] ). The present approach shows that the GP equation 
with an arbitrary potential h(p) can be derived from the 
KG equation with a self-interaction potential V (p). More 
precisely, Eq. (C14| shows that the potential h(p) in the 
GP equation is equal to the derivative of the potential 
V(p) in the KG equation. Reciprocally, V{p) is a prim¬ 
itive of h{p). The primitive of h(p) played some role in 
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our previous studies [55], and was noted H(p). When the 
GP equation is derived from the KG equation, we have 

H(p) = V{p). (C15) 


The KG equation can be written in terms of h as 


1 d 2 4> 7 

We also have 


h 2 


1 + 


2 $ 


vm 2 ) = h 


, m 7 


m 


n 2 


h 2 


I# 


b = 0. 

(C16) 

(C17) 


3. The Madelung transformation 


Using the Madelung 90 transformation 

%l> = Jpe iSlh , u = — VS, (C18) 

m 

where p{ r, t) is the density, S(r, t) is an action, and u(r, t) 
is interpreted as an irrotational velocity field, we can 
rewrite the GP equation ( |C13[ ) in the form of hydrody¬ 
namic equations (see, e.g.. [89] l: 

^+V-(pu)=0, (C19) 


+ (u • V)u = —-VP - V4> 
at p 

where 

v 2 m y/p 


-VQ, (C20) 
m 


(C21) 


is the Bohrn quantum potential and P is the pressure. 
The pressure is given by a barotropic equation of state 
P = P(p) determined by the potential h(p) according to 


h'(p) = 


P'(P) 


(C22) 


Since \7h = (l/p)VP, the potential h in the GP equa¬ 
tion has the interpretation of an enthalpy in the quantum 
Euler equation (C20). We note that the pressure is ex¬ 
plicitly given by 


P(p) = ph(p) - H(p) = pV'(p) - V(p). (C23) 


4. Energy functional 


Th e barotropic quantum Euler equations (C19) and 
(C20) conserve the mass and the free energy functional 27 


F = E + U, 


(C24) 


2 ' As a result, a minimum of free energy F at fixed mass M is a 
stable stationary state of the GP and quantum Euler equations. 


where E = 0 C + ©q + W is the sum of the classical 
kinetic energy O c , the quantum kinetic energy Oq and 
the potential energy W. and U is the internal energy. 
They are defined by (see, e.g., [55]): 


Qr = 


p Y dr ’ 


Bq = - 
m 


pQ dr, (C25) 


W = J dr, U = Jpf dp 1 dr. (C26) 

The internal energy can be written as 

U = J u(p) dr, 


(C27) 


where 


l (p) = P j 


P(P') 


dp' 


(C28) 


is the internal energy density. 

The internal energy density u(p) defined by Eq. (C28) 
corresponds precisely to the term that appears in the en¬ 
ergy density of Eq. (A15) in addition t o the rest-mass 

depending 


IIH 


density pc 2 . As we have indicated in Sec. 
on the pressure law P(p), this term mimics a “new fluid” 
that adds to “dark matter”. This new fluid may be an ex¬ 
otic constituent (e.g. a stiff fluid) appearing in the early 
universe [44]. It may also represent the “dark energy” 
in the late universe (see Sec. IIHi. Using the present 
formalism, we can make a connection between this new 
fluid and the potential that appears in the GP and KG 


equations. Integrating Eq. (C28) by parts and using Eq. 


(C22), the internal energy density may be rewritten as 

u{p) = ph(p) - P{p). 

(C29) 

According to Eqs. (|C15), (|C23) and (C29), 

we obtain 

u{p) = H(p) = V{p), 

(C30) 

or, equivalently, 


(cT 

II 

(5) 

II 

V 

jT 

II 

(C31) 


Therefore, the internal energy density u(p), which mimics 
a “new fluid” in Eq. (A151, is equal to the potential V(p) 
appearing in the KG equation. 28 

Remark: We note that H(p) and V(p) are defined up 
to a term of the form ap + b. Therefore, we can adapt 
the coefficients a and b in order to have the simplest 
expressions of H and V. We will use this prescription in 
the following section. 


It is determined by the Euler-Lagrange equation 8F — fi5M = 0, 
where fi is a Lagrange multiplier (chemical potential), leading to 
+ mh(p) + Q = mp, which is equivalent to the condition of 
hydrostatic equilibrium VP + + (p/m)VQ = 0 |89| . 

28 Actually, this equivalence is valid only in the nonrelativistic limit 
c —¥ +oo. The relativistic regime is more complicated to investi¬ 
gate (since p is not the rest-mass energy) and will be considered 
specifically in a future communication. 
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5. Particular examples 

The internal energy density u(p) is determined by the 
equation of state P(p). As we have seen, this equation of 
state can be obtained from a field theory based on the KG 
or GP equation. In this section, we consider particular 
examples of equations of state. 


a. Isothermal equation of state 

We consider the isothermal equation of state [Z2: 

k B T 


P = p- 


(C32) 


In the present context, this equation of state is expected 
to arise from a self-interaction potential, not from ther¬ 
mal motion. As a result, the temperature has to be re¬ 
garded as an effective temperature which can be either 
positive or negative. The internal energy is 


k B T . 

u = - pinp. 


The potential in the GP equation is 

k B T 


, k B T 

ri = -m p, 

m 


H = 


m 


P hi P- 


(C33) 


(C34) 


The GP equation takes the form 


dib k 2 

ih-—- = - Aib + m&il) + 2k B Tln \ib\ib. 

at 2m 


The free energy can be written as F = E—TS where E is 
the energy, T is the temperature and S is the Boltzmann 
entropy 


S = —k B / — In — dr. 
J m m 


The potential in the KG equation (Cl I is 


K(H) = 2 ^^|</>| 2 lnH, 


(C36) 


(C37) 


so the KG equation writes 

1 9 2 </> m 2 c 2 

ct> + ~h 2 ~ 


2 $ 
1 4 — 2 


4 mk B T f m\4>\ 


= 0 . 


(C38) 


where the polytropic constant K may be positive or neg¬ 
ative (for the same reason as before), and the polytropic 
index 7 is arbitrary. The internal energy is 


K 

7-1 


p 1 = 


P 


7-1 


(C40) 


Since P = ( 7 — l)u, this component in Eq. (32) can mimic 
dark energy when 7 —» 0 (see Sec. IIH). For 7 = 0, the 


pressure is constant (P = K) and we recover the equation 
of state P = — u of dark energy. The potential in the GP 
equation is 


h = 


TV- 1 , 


H = 


7-1 

The GP equation takes the form 


K 

7 - 1 / 


(C41) 


= — AV' + + K{n + (C42) 

(J L Z / / L 

The free energy can be written as F = E — KS„, where 
E is the energy, K is the polytropic temperature and S~ f 
is the Tsallis entropy [53] : 


5 

1 

7-1 


j(P' ~ P) dr. (C43) 


The potential in the KG equation (Cl I is 

K / TO \ 




27 


(C44) 


(C35) so the KG equation writes 


c 2 dt 2 ^ 


2 2 
m A c z 


h 2 


, 2 $ 

1 -I-w 


+ 


2m 2 K r y f m 


2(7-1) 


(7 — 1 )H 2 ly h 
The usual GP equation writes m-- 

4-7T dell 2 


= 0 . 


(C45) 


= 

dt 2m 


■ 


IV’lV (C46) 


It describes a BEC at T = 0 with a q uadra tic self¬ 
interaction {h oc |^| 2 ). The GP equation (C46) can be 
derived from the mean held Schrodinger equation with 
a pair contact potential (see Sec. II.A. of [55]'). In that 
case, a s represents the scattering length of the bosons. It 
is positive when the self-interaction is repulsive and nega¬ 
tive when the self-interaction is attractive. The potential 
in the GP equation (C46) is 


b. Polytropic equation of state 


47 Ta s H 2 

}l = P. 


H = 


2ira*h 2 


~P 


(C47) 


We consider the polytropic equation of state [77]: 


P = Kp\ 7 = 1+-, 


(C39) 


The pressure is 


P = 


2ira s h 2 


(C48) 
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It corresponds to a polytrope of index n = 1 (i.e. 7 = 2) 
and polytropic constant K = 27ra a 7i 1 2 /m 3 4 5 6 7 . The internal 
energy is 


u = Kp 2 = P. 


(C49) 


Since P = u, this component in Eq. (32) can mimic stiff 


matter [44] . Historically, the stiff equation of state was 
introduced by Zel’dovich issuszi in the context of baryon 
stars in which the baryons interact through a vector me¬ 
son field. The free energy can be written as F = E—KS 2 
where E is the energy, K is the polytropic temperature 
and S 2 is the Tsallis entropy of index 7 = 2 : 


S 2 = - P 2 dr. 


The potential in the KG equation (Cl) is 




K 2 


so the KG equation writes 


A, 

c 2 dt 2 ^ 


2 2 
me 


h 2 
87 t a s m 
h 2 


2 $ 
1 4 - 7T 


m=o. 


(C50) 


(C51) 


(C52) 


Since P = —u — A , this component in Eq. (321 can mimic 


dark energy. The potential in the GP equation is 
A 


h = — 


H = -AliL(p/p*) - A. (C55) 


The GP equation takes the form 

dib h 2 1 

ift— = — -—Atp + m&i/j — Am—^ip. (C56) 


dt 


2 m 


h/f 


Eq. (C56) can be viewed as a GP equation with an in¬ 
verted quadratic potential, i.e. with the exponent —2 


instead of +2 in the usual GP equation (C46). To our 


knowledge, this equation has not been introduced before. 
This equation can be obtained as the limiting form of Eq. 
(C42) when 7 —► 0 (n —> —1) and K —>• 00 with A = A '7 


finite. However, it cannot be obtained from a polytropic 
equation of state with 7 = 0 (which corresponds to a 
constant pressure P = K) because a constant pressure 
yields h = 0 (or h = cst.) in the GP equation, which is 
different from Eq. ( C56| ). This shows the “regularizing” 
property of the logotropic equation of state when 7 —> 0 , 
as in the case of the Lane-Emden equation (see Sec. 0- 
The free energy can be written as Fl = E — ASl where 
E is the energy, A is the logotropic temperature and Sl 
is the log-entropy [55] : 


c. Logotropic equation of state 


Sl= lap dr. 


(C57) 


We consider the logotropic equation of state 
P = Aln(p/p*), 


(C53) 


where the logotropic constant A may be positive or neg¬ 
ative. The logotropic equation of state (C53) can be 


viewed as the limiting form of the polytropic equation 
of state (C39) when 7 —>• 0 (n —> —1) and K —► 00 in 

The internal 


such a way that A = K'y is finite 
energy is 
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u = — Ahi(p/p*) — A = —P — A. 


(C54) 


The potential in the KG equation ( |Cl| ) is 
V(\cP\) =-2Aln\<P\, 
so the KG equation writes 


(C58) 


1 d 2 (p _ 

c 2 dt 2 * 


H 2 


2 $ , 

1 + -s- <£- 


2 A 


= 0 . 


(C59) 
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